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[1] Learning Modelling with DERIVE , Townend, Poutney 
Prentice Hall, 244 pages, ISBN 0-13-190521-X 

[2] ACTAS de las Jornadas sobre la ensehanza de matematicas con Derive 

You can find some information in the User Forum, page 5. 

[3] Der TI-92 im Unterricht, Klaus Aspetsberger u. Franz Schloglhofer 

[4] Mathematik erleben mit dem TI-92, Gunter Schmidt 

Both publications are available from Texas Instruments Deutschland, 85356 Freising (Tel:08 161 804984) 


Exchange for DERIVE Teaching 


Borse fur DERIVE Unterrichts- 

materials in the DNL 


materialien im DNL 

The wheel has not to be invented twice. 


Das Rad muss nicht zweimal erfunden werden. 


I can offer: 

Binomial Theorem, GCD & LCM, System of coordinates, Modeling Word problems with DERIVE 
(all in English and German), SET.EXE, MENGE.EXE. Functions -Domain and Range from Tom 
Drummond, Glasgow. I have produced a paper "Einfuhrung in die Matrizenrechnung mit DERIVE ' -in 
German. If you are interested I would send you the paper on a diskette in MS-Word6-format. 

In the last time I received many interesting contributions to be published in the next DNLs. They all 
were accompanied by friendly letters. I'll reprint some sentences to give an idea about their contents: 

David Halprin, North Balwyn, Australia: ... I have opened a mathematical 'can of worms', too much for one 
person to develop to its full potential, so I invite fellow members of DUG or readers of the Newsletter to write to 
me (about the "Cesaro Glove-Osculant"). 

Maria Koth, Vienna, Austria: . . . Motiviert durch unser Gesprach auf der Lehrerfortbildungstagung habe ich 
jetzt meinen Artikel „Computergrafik mit DERIVE 44 ("Computer Graphics with DERIVE") iiberarbeitet und 
mochte Ihnen die verbesserte Version zusenden . . . (about generators of fractals, Sierpinski triangles a.o.) 

Prof. M.J.Fernandez Guitierez, Oviedo, Spain: I now enclose copies of our manuscript "Solving third-order 
linear differential equations with constant coefficients", that I believe fits the requirements for papers in the 
DERIVE NEWSLETTER. So I would like you to consider it for publication . . . 

G P Speck, Wanganui, New Zealand: The DNL readers may find the catalogue of function tables presented 
with attendant MTH files a significant time saver as I have on several occasions. 

Peter Mitic, Medstead, England: I enclose three articles which you might look at. They are an extension of the 
talk I gave at the Plymouth conference in 1994, "Exploiting new Features in DERIVE 3", "The Normal Distribu- 
tion: two Proofs and a Simulation", "Probability Distributions: Proof and Computations" . . . 

Many thanks for your long and contentful letter. Parts of it will be the introduction to your contribution. 

I would like to receive your ’short notes' on the use of probability generating functions. I think that 
probability theory is a very important part of mathematics and we could enforce it in the DNL. So I 
published one teaching unit in this DNL which I wrote some time ago for use in my classes. Once 
more many thanks and much luck in your new career. Josef. 

Prof. Neil Stahl, Menasha, Wisconsin, USA: . . . For years I felt computer graphics would be of benefit in dem- 
onstrating important concepts in calculus, such as tangent lines and tangent planes as well as direction fields. I 
am pleased at the opportunity you are offering me to distribute this work more widely. 

Nurit Zehavi, Rehovot, Israel: I hereby submit a "DERIVE tip" to the DUG newsletter. I hope that you and 
other users will find it helpful. 


Many thanks to all of you. I think we all can be proud of the intemationality of our group and of its en- 
thusiasm , of course. We cooperate from country to country and from contintent to continent. Let's pro- 
ceed on our way! 
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Liebe DUG Mitglieder, 

Ich habe versucht, in diesem DNL einige alte 
„Schulden“ zu begleichen, d.h. Beitrage aufzuneh- 
men, die schon lange auf meinem „ Stack “ liegen. 
Anstelle der "TITBITS" finden Sie dieses Mai einige 
Anmerkungen und Erganzungen zu den letzten Tit- 
bits. Ich finde den e-mail- Austausch besonders in- 
ter essant, weil uns Albert Rich und Johann Wie- 
senbauer damit einigen Einblick in das Innenleben 
von DERIVE ermoglichen. 

Ich danke Alfonso J. Poblacion fur die 
Anregung zum AC DC-Teil und auch 
dafur, dass er gleich den ersten Beitrag 
dafur geleistet hat. Ich mochte Sie gerne 
einladen, sich fallweise an der Gestaltung 
dieser - vielleicht standigen - Kolumne 
des DNL zu beteiligen. 

Fur die TI-92 Benutzer und solche, die es 
werden wollen, ist sicher interessant, dass 
ein von TI unterstiltztes Fortbildungsprogramm I s 
(Teachers Teaching with Technology), das in den 
USA unter der Leitung von Bert Waits und Frank 
Demana schon etabliert ist ,auch in Europa als 
1 6 -Europa eingefuhrt werden soil. Ich habe die Ge- 
legenheit, im Juni in Columbus, OHIO, an einer 
I s -Sommerschule teilzunehmen und kann im nach- 
sten DNL dariiber berichten. 

Nun zwei Hinweise: Ihre DUG Mitgliedsnummer 
finden Sie auf dem Adressenaufkleber, und mit glei- 
cher Post sende ich die MTH-files zu diesem DNL 
nach Hawaii. Albert Rich hat zugesagt, diese nach 
Moglichkeit in seine web-Seite aufzunehmen. Viel 
Gluck beim „Downloaden “! 

(http fwww. derive, com). 

Wir freuen uns schon auf Bonn. Sie konnen als Vor- 
geschmack drei mit DERIVE berechnete Raumfigu- 
ren sehen. die ich mit einem Hilfsprogramm ins 
ACROSPIN-Format konvertiert habe. Damit kon- 
nen diese Figuren nun auch animiert werden. 

Mit den besten Grilfien, Ihr 
Josef Bohm 



Conic sections 


Dear DUG Members, 

I have tried to settle some old debts in this 
DNL, i. e. to publish contributions which have 
been waiting on my "stack" since long. Instead 
of the TITBITS you will find some completions 
to earlier articles. I find the e-mail exchange 
especially interesting because Albert Rich and 
Johann Wiesenbauer give some insight into 
inside DERIVE. 

My special thanks go to to Alfonso J. Pobla- 
cion not only for his idea to create AC 
DC but also for submitting the initial 
paper. I would like to encourage you to 
contribute to this - maybe permanent - 
new column. 

Certainly it is interesting for being and 
ongoing TI-92- Users, that T 3 
(Teachers Teaching with Technology) - 
a profess-ional development program, 
spons-ored by TI, which was founded 
by Bert Waits and Frank Demana - will be es- 
tablished as T 3 -Europe in Europe. 

Fortunately I have the opportunity to attend the 
next f summer school at Columbus, 0/-//0 in 
this June. I will give a report about this event in 
the next DNL. 

Two notes: You can find your DUG member- 
ship number on your address label, and I am 
mailing all the MTH-files belonging to this DNL 
issue to Hawaii. Albert Rich has promised to 
put them on SWHH's website. Much luck with 
downloading. (http://www. derive, com) 

We are looking forward to the Bonn Confer- 
ence. As a foretaste you can see three DE- 
RIVE calculated 3D-objects which I converted 
into the ACROSPIN format using a self made 
program. It enables animating those figures 
using different colours and layers. 

With my best regards 
Josef 



Elliptic torus with a torus knot line and its tube 
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The DER1VE-NE WSLE TIER is the Bulle- 
tin of the DERIVE User Group. It is pub- 
lished at least four times a year with a con- 
tents of 30 pages minimum. The goals of 
the DNL are to enable the exchange of ex- 
periences made with DERIVE as well as to 
create a group to discuss the possibilities 
of new methodical and didactical manners 
in teaching mathematics. 

We include a section dealing with the use 
of the TI-92. 


Editor: Mag. Josef Bohm 
A-3042 Wiirmla 
D'Lust 1 
Austria 

Phone: 43-(0)660 40 70 480 
email: nojo.boehm@pgv.at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English ar- 
ticles are very welcome nonetheless. Your 
contributions will be edited but not as- 
sessed. By submitting articles the author 
gives his consent for reprinting it in the 
DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE Newsletter will be. 


Preview: (Contributions for the next issues): 

Graphic Integration, Linear Programming, Bohm, A 

LOGO in DERIVE, Lechner, A 

3D Geometry, Reichel, A 

Parallel- and Central Projection, Bohm, A 

Algebra at A-Level, Goldstein, UK 

Tilgung fremderregter Schwingungen, Klingen, GER 

Utility for Complex Dynamic Systems, Lechner, A 

Some notes on DERIVE 2.6 functions and limits, Speck, NZL 

Linear Mappings and Computer Graphics, Kiimmel, GER 

Solving Word Problems with DERIVE, Bohm, A 

DERIVE and ACROSPIN, Schorn & Bohm, A/GER 

Visualizing a Special Line in the 3D-Space, Zehavi, ISR 

Line Searching with DERIVE, Collie, UK 

The TI-92 Section, Bert Waits, Bernhard Kutzler, Frank Demana 
and 

Setif, FRA; Vermeylen, BEL; Leinbach, USA, Aue, GER; Halprin, AUS; 
Weth, GER; Wiesenbauer, A; Keunecke, GER; Weller, GER; ... 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
Herstellung: Selbstverlag 
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Pierre A. Arnoldi, Vermont Switzerland 

Derive for Windows will no doubt be of great interest to me. I am running an IBM compatible machine with 
EGA/VGA screen. My printer is a CANON BJC-70 (black and white and four colors printing head). I was not 
yet able to connect some fake-printer to the machine for my graphs to be printed in colors. Can you tell me 
which printer I should choose among all printers registered in Windows, which works the same way my Canon 
works (ink bubble system)? How am I supposed to select the DERIVE printing options then? 

... You should note that I am not a mathematician at all. I happened to study some algebra when I was young as 
well as geometry and other related branches of such sciences. I was more fascinated than really interested and I 
am not going to be another Einstein . . . 

DNL: I asked Micheal Petz from SWHH about your problem. Michael has a Canon BJ for his own (only black 
& white) and he is running his machine in Epson mode. He does not face any problem. I have tried with an HP- 
Deskjet Colour printer and had no problems, too. You can find the appropriate DeskJet C setting in the Transfer 
Print Options. Much luck for producing nice colour prints with DERIVE. 

John Berry, Plymouth, UK 

Dear Josef, HELP!! 

You are the DERIVE wizz, can you solve my problem below? 

This set of steps produces a part of the bifurcation diagram for the quadratic map x - c. I want to automate it for 
lots of values of c, plotting the answers. Any idea how I can do it please? 


ORBIT (v 

#1: 

i g, x, n, c) 

#2: 0RB_DIAG(g, x, a, n, 

#3: 1 

F(c, n) 

:z ORBJDIAGCx 

#4: 1 

F(1.B, 

53) 


' 1.5 

0.615312359 

#5 : 

1.5 

-1.1213907 


. 1.5 

-0. 2424S2S96 


:= VECTOR (|Y , v j, r, El, 
c) := ORBIT(ITERATES(g , x 


2 

- c , x , 0, n , c) 


n) 


n) , g, x, 


n , 


O 


DNL: I hope that the “wizz” can help. If not then please don 7 blame my wizzardry. Proceed with: 

MQRE_F(n, start, end, inc) := VECTOR (F(c, n) , c, start, end, Inc) 
M0RE_F(E3 , 0, 2, 0.01) 


■2 

1.5 

■1 

0.5 


-0.5 

-1 

-1.5 


■ . . . . . . . . . . u . , 

4 _ 0.5 0 . 6 0.7 0.,>"tr5"' 1 1.1 1.2 1.3 "l;4 a 1,5 1 . 6 . 1 . 7 ' 1 . S- 1 . 9 2 

■■ , ■ '■ . 


-2 
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You will find a Tl-implementation of this algorithm in the TI-92 Comer. Josef 

Helmut asked for tricks. Here you are: The QUOTIENT-f unction works with real coefficients, so let's 
disguise the complex i. I replace the complex #i by a real variable i, perform the divisions and then 
undo the former substitution. Let's look if I can outwit DERIVE 6? 

Let's produce a quotient function which works with complex numbers, too: 


#30: 

com| 

#31: 

com| 

#32: 

p2 : 

#33: 

p2 : 

#34: 

p3 : 

#35: 

P3 

#36: 

p4 

#37: 

p4 


= z + 5. 69384- z - 1,5922 + i- (1,62108- z + 9,61268) 

= compt_quot(p3 , z - z3) 

= z + 5,95897 + 0,007596 - L 

As you can see, it works fine. 

Another possibility is to Expand Trivial p(z) and then substitute all remaining fractions by 0: 

p 


#38: 


z - zl 

Expand trivial content 

-7 li li 

2.68244-10 -(9. 91872-10 ■ z - 6 . 59019-10 ) 3 2 

#39: + z + 5.4783- z - 0. 395571- z + 14.7164 + 

11 2 11 10 

2,5-10 -z -2.3915-10 - z + 6 .11532-10 

-7 11 11 

5.85026-10 -(6.72502-10 ■ z - 2 .64415-10 ) 2 

+ 0.125866- z + 0. 749733- z + 0.308808 


11 2 11 10 
2.5-10 -z -2.3915-10 ■ z + 6 .11532-10 


3 2 2 

#40: p2 := z + 5. 4783- z - 0. 395571- z + 14.7164 + L- (0.125866- z + 0. 749733- z + 0.308808) 


and so on ... 


Alfonso J. Poblacion Saez, Valladolid, Spain 

Dear Josef, 

I am glad to write to you again and absolutely fascinated about the three last Newsletters, 
congratulations. 

In last June a work meeting group about teaching mathematics with DERIVE took place in 
Santander, Spain, for which I was a co-organizer. So, I send you the proceedings of this Spanish 
meeting. I realize that it is a little late but it was impossible for us to have them ready earlier. The 
company which sponsored the meeting (Technical Research from Barcelona) promised to edit the 
proceeding, and after waiting and waiting (as patiently and innocently as you can imagine) we finally 
decided to do the editing by ourselves. It is not a marvellous edition as you can see but it is not too 
bad at last. 

Unfortunately they are in Spanish except the last dossier of the addendum about the TI-92 by 
Bernhard Kutzler who introduced it to us. We were some of the first lucky people who handled the 
machine. I also include two pages, one with a translated summary of the aims of the meeting (this is 
- more or less - the foreword on pages 4 to 6 of the proceedings) and the English description of the 
following contents (pages 1 to 3). 
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I hope you or some of the DUG members will find them useful. If you or someone is interested in any part and 
has problems with the language, I can translate it into English for you/him/her. 

My e-mail is alfonso@gauss.mat.eup.uva.es. I prefer, if it is possible this way of communication. Thanks for the 
patience of reading this letter. 

Hope to hearing from you soon. Yours faithfully A.J.P.S. 

FOREWORD (not exactly; I summarize and include some other information) 

In 1992 a group of Spanish mathematics teachers at University level began to interchange information 
and experiences about the use of DERIVE, mainly in the first courses of technical careers. In Septem- 
ber of 1 993, this group set up a Spanish DERIVE Users' Group with the aim of organizing and streng- 
thening activities about teaching with DERIVE. As we also think that the wheel has not to be invented 
twice, the Group sent to each member all the information got and required. The Group has no eco- 
nomic resources, so we established an itinerant center of operations per year to share expenses and 
time. Each year it sends information at least four times. The first center was Madrid (U.P.M.) and the 
chief responsible was Rafael Minano. The second one was the Department of Matematica Aplicada a 
la Tecnica (Applied Maths to the technic is the approx, translation) at the University of Valladolid and 
Alfonso J. Poblacion and Carlos Marijuan were the persons in charge. Today the center is in Valencia 
and Jose Luis Llorens is the main responsible. 

After the Spanish presence in the First International DERIVE conference in Plymouth, our Group 
thought in the convenience of organizing a work meeting where a few number of selected participants 
could discuss deeply in the incidence of DERIVE in teaching. On page 4 of the proceedings you can 
find the Organizing Committee that pointed out the discussion themes and chose the experienced par- 
ticipants. They were divided into three groups: 

Group I: Algebra, Geometry and Discrete Mathematics 
Group II: Calculus and Numerical Analysis 
Group III: Mathematics at Secondary School 

Each group worked independently with different moderators per work session (you can see them on 
pages 5 and 6). 

Short overview of the Contents: 

1. Works of Group I: Scheme of work, Introduction and Experiences, Methodology (General ideas and ideas 

how to use DERIVE day by day), Evaluation, Incidence on learning and in curricula, Benefits and disad- 
vantages about using DERIVE. 

2. Works of Group II: Scheme of work, Experiences, Objectives and Methodology, Incidence on learning, Influ- 

ence in curricula, Evaluation with DERIVE. 

3. Results of a survey answered by Group I and II. 

4. Works of Group III: Scheme, Description about the particular situation in these studies, Methodological con- 

siderations, Some conclusions, Results of a survey among the participants and among other secondary 
school teachers. 

5. Comparison between DERIVE and other systems 

6. Proposals for future DERIVE versions 

7. Conference: Calculus Discretization. Picard's Method with DERIVE by Miguel de Guzman 

8. Addendum: Dossier Group I, Dossier Group II, Dossier III: TI-92 

DNL: Many thanks , Alfonso for your paper (190 pages). Many thanks also for your generous offer to 
translate any part. I hope you will find some interested DUG members who will contact you. (But I also 
hope that there will not be too many!!!). I would like to say that especially the dossiers are very useful, 
because they consist of some materials provided by each participant on which you based your dis- 
cussions. Once more many thanks and good luck for the Spanish DUG. I hope that we will be able to 
enforce our cooperation. I have the pleasure to confirm that we have a remarkable increase in Span- 
ish members in the DERIVE User Group. Ole y much as gracias.l 

One more thank for your splendid idea of AC DC. You will find it at another place in this issue. 


Find more USER FORUM activities questions on the last page. 
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IMP Spider and SPIRAL GRAPH 

Edward Sawada, Mililani, Hawaii 

It is amazing where you can find a math problem. At an interactive Mathematics Program 
(IMP) workshop in San Francisco, I noticed the logo on the IMP handouts was the spiraling 
square. This logo is illustrated in the book The Joy of Mathematics on page 228, Spider and 
Spirals. I took it upon myself to see if I could duplicate a similar logo using DERIVE. The 
following is my effort to show the mathematics used in trying to duplicate a similar spiral. 

To make maximum use of the 1 to 1 scale for DERIVE 's Plot Window, I chose to use the 
range of values for x and y as [-3,3], and I decreased each side of the square by 0.1. A minor 
problem was to figure the approximation for the angle of rotation for each square. I resolved it 
in the following manner: 

O.lx 

SIN(a) = 0. 1/0.9 
a = 6.37937° 

The next problem I faced was to find a way to reduce and rotate the matrix of points that 
made up the initial square. 

#1: CaseMode := Sensitive C-3,3)^ C3,3) 

'3 3 

3 -3 

#2 : square := -3 -3 


From analytical geometry we can verify the rotation formula involving the point (p,q) ro- 
tated by cf as (p cos(a) - q sin(a), q cos(a) + p sin(a)). 

ROTATION OF A POINT 

p q 

#3: COS(p) = , SIN(p) = 

r r . 

#4: x — r-COSCa + p) 

#5: Trigonometry := Expand 

#6: x = r-CQS(a) -COS(p) - r- SIN( h) ■ SIN(p) 

P q 

#7: x — r-COS(a) r-SIN(a) 

r r 

#8: x = p-CQS(a) - q-SIN(a) 

#9: y = r- SIN (a + p) 

#10: y = r-COS(a)-SINCP) + r- SIN (a) -COS(p) 

r-C0S(a)-q r-SIN(«)-p 

#11: y = + 

r r 

#12: y = q-COS(a) + p.SIN(a) 
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From the logo, I could see that each square was reduced and rotated. The challenge there- 
fore was to write a program that does just that. The ITERA TEN-command allowed me to use 
the recursion technique to reduce and rotate as many times as I wanted. The points to be ro- 
tated were identified by the ELEMENT-command for matrix operation and rotated with the 
rotation formula, r denoted the reduction ratio, 8 stands for the degree of rotation, and n will 
be the number of rotations desired. 

(Comment from 2009: Instead of ELEMENT we use now the SUB-functionality.) 


fT r-a -C0S(5) - r- a 
1,1 1 

r-a -C0S(5) - r-a 
2,1 2 

r-a -C0S(5) - r-a 
i'rnpspi'der_web(r , 5, n) := ITERATES 3,1 3 


r-a -C0S(5) - r-a 
4,1 4 

r-a -C0S(5) - r-a 

IL S,1 3 


SIN (5) 

r- a 

■C0S(5) 

+ 

r-a ■ SIN (6) ' 



1, 

,2 


1,1 


SIN (5) 

r- a 

■CQS(5) 

+ 

r-a ■ SIN (6) 



2, 

,2 


2,1 


SIN (6) 

r- a 

■CQS(5) 

+ 

r-a ■ SIN (6) 



3, 

,2 


3,1 

, a, square, n 

SIN (6) 

r- a 

■CQS(5) 

+ 

r-a ■ SIN (6) 



4, 

,2 


4,1 


SIN (6) 

r- a 

■CQS(5) 

+ 

r-a ■ SIN (6) 



3, 

,2 


3,1 

J 


This is a more elegant version: 


fmpsp_web(r, 5, n) := ITERATE S(VECT0R(T r- a -C0S(5) - r-a ■ SIN (5) , r-a -C0S(5) 

L n ',1 i ,2 i ',2 


r-a -SIN(S)"|, i, B) , a, square, n) 

i,l J 


If you prefer working with matrices then do it the following way: 


nota(S) 


‘ C0S(S) 
. SIN(S) 


- SIN(S) ' 
COSCS) _ 


l'rnpjiatCr, S, n) := ITERATES(r- a- nota(S) , a, square, n) 


IMP SPIDER WEB 



CORNER SPIRALS 



impspider_web(0.9, 6,37937°, 30) C0PR0]ECTI0N(inipspider_web(0,9 , 6,37937°, 30)) 

nmpsp_web(0.9 , 6,37937°, 30) C0PR0]ECTI0NCin,psp_web(0,9 , 6,37937°, 30)) 


Lesson Learned 

1. Comer spirals created by 'transposing' matrix (using the COPROJECTION function). 

2. This program will help me in setting up future programs where I need to reduce and ro- 
tate a matrix of points or parametric equations. 


plO 
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{Two Weeks Later} 

I am not satisfied by my creation of the IMP WEB, for I know I could have done a better 
job if I had spent more time on it. Now I have the time, so I here present another method to 
create the IMP WEB and a graph I will call the Misguided Missiles. The technique used in 
both graphs are very similar. The method I used comes from a lesson I once presented on vec- 
tors. The lesson was how to divide a segment into ratio of a : b using vectors, the end points 
of the segments are known. 



A = x 1 +i-y l , B = x 2 +i- y 2 , C = B-A, C, = C 

a + b 

C = x 2 +i-y 2 -(x l +i-y l ) = x 2 -x l +i(y 2 -y l ) 

a a ( , u ax 2 -ax l +i(ay 2 -ay l ) 

c ' = ^- b c =^-x, + '(yi-y,))= — b 

A i n _(x l +i-y l )(a + b) i ax 2 -ax l +i(ay 2 -ay l ) _bx l +ax 2 i i(by l +ay 2 ) 

L) — A + L | — I — I 

a+b a+b a+b a+b 


bx, +aXr, 


by,+ay 2 


The x and y vector component will be used to increment sides of a square or an equilat- 
eral triangle by a ratio a : b. (x, y) is the new coordinate after the division. The ITERATES 
command was used for subsequent reduction in the sides of the polygon. By increasing the 
value of n, I can increase the spiraling effect. By adjusting the ratio a : b I can vary the direc- 
tion of the spiral. 
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IMP WEB 



IMP WEB 


MISGUIDED MISSILE 



MISGUIDED MISSILES 





i mp_web2 (1,1, square , 10) 


imp_web2(l, -20, square, 10) 


MISGUIDED MISSILES 

o ^ " 

2-jl -3 

tn 

- 2-^3 


L o 3 J 
imp_web2(40, 3, tri, 20) 



0 

3 " 


0 

-3 " 


" - ft 

0 ' 


" ft 

0 ' 



ft 

0 


ft 

0 


0 

-3 


2' ft 

-3 


tnl := 

- ft 

0 

, tri2 := 

- ft 

0 

, tri3 := 

- 2-ft 

-3 

, tri4 := 

0 

-3 



0 

3 . 


0 

-3 _ 


- - ft 

0 . 


- ft 

0 . 



[ imp_web2(40 , 3, tnl, 20), -imp_web2(40 , 3, tn2, 20), imp_web2(40, 3, tn3, 20), imp_web2(40, 3, tn4, 20)] 

Mathematics is the science of patterns. The patterns maybe number patterns or geometric 
patterns. Here, in the graphs I have created, I used number patterns to get the geometric pat- 


terns. 
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Because of programs like “ DERIVE ”, I firmly believe that the future of mathematics 
will be in mathematical analysis and the programming of mathematics. 




' 0 

0 ' 


" -3 

0 " 


' 0 

0 ' 


" 0 

0 " 



3 

0 


0 

0 


-3 

0 


0 

-3 


sql := 

3 

3 

. *q2 := 

0 

3 

. := 

-3 

-3 

. sq4 := 

3 

-3 



0 

3 


-3 

3 


0 

-3 


3 

0 



. 0 

0 . 


_ -3 

0 _ 


. 0 

0 . 


_ 0 

0 _ 



[ imp_web2(4 , 1, sql, 10), nmp_web2(l, 4, sq2, 10), imp_web2(4 , 1, sq3, 10), nmp_web2(l, 4, sq4, 10)] 

Here you can find a generalization of the problem to create other patterns. Josef 




Rotating and Shrinking a Polygon 
#35: rot (a, q) := q-rota(a) 



r 

f 2 ■ 7T ■ k 1 


#36: polyCr, n) := VECTOR 

rot 

, [r, 0] 

, k, 0, n 



^ n j 

j 



#38: pol_rot(r, n, cp, j) := ITERATE5(s(n , tp)' ro ' t (tpi p) , p, poly(r, n) , j) 



pol_rot(3, S, 10° , 20) 


C0PR0JECTI0N(pol_rot(3 , S, 10° , 20)) 
C0PR0JECTI0N(pol_rot(3 , S, - 10° , 10)) 
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Ebene Algebraische und 
Transzendente Kurven (8) 

Thomas Weth, Wurzburg, Germany 

Didaktischer Nachtrag zur Konchide des Nikomedes (Folge 5) 


1. Einleitung 

Im herkommlichen Geometrieunterricht werden 
mittlerweile im Allgemeinen nur noch geradli- 
nig begrenzte Figuren und der Kreis behandelt. 
Die Abbildungen, die im Unterricht thematisiert 
werden, sind ausschlieBlich Ahnlichkeitsabbil- 
dungen, die Geraden auf Geraden, und Kreise 
auf Kreise abbilden. In neuester Zeit werden 
unter dem Einfluss von Computerprogrammen 
wie Cabri-Geometre Versuche untemommen 
(z.B. Werge/Bock), „wenigstens“ die Inversion 
am Kreis in den Unterricht einzubeziehen. Aber 
auch die Inversion ist eine zykel-treue Abbil- 
dung, die die Menge der Geraden und Kreise in 
sich selbst abbildet: „echte“ Kurven werden 
beim Abbilden von Kreisen und Geraden nicht 
erzeugt. Neben der Formenarmut des Geomet- 
rieunterrichts ist eine zweite Schwachstelle die 
Trennung zwischen Algebra und Geometrie. 
Alleine schon die allgemeine physikalische 
Trennung von Algebra- und Geometriebuchern 
und -heften suggeriert dem Schuler nachweis- 
lich (vgl. Weth 1993), dass Algebra und Geo- 
metrie zwei disjunkte mathematische Diszipli- 
nen sind. An dieser Einstellung hat auch die 
Klein'sche Reform des Mathematikunterrichts 
nur wenig geadert. Klein hatte versucht, dem 
Schuler iiber die Behandlung von Funktionen 
die „mathematische Wissenschaft als ein groBes 
zusammenhangendes Ganzes 66 zu vermitteln, da 
iiber die Darstellung von Funktionsgraphen der 
geometrische und iiber die Behandlung und 
Diskussion von Funktionstermen und - 
gleichungen der algebra-isch/analytische As- 
pekt gleichberechtigt behan-delt werden konn- 
ten. 

Die beschriebene Trennung von Algebra und 
Geometrie zieht sich mittlerweile auch bis in 
die gymnasiale Oberstufe, wo die ehemals ana- 
lytische Geometrie zu einer sterilen „Linearen 
Algebra 46 degeneriert ist. Immer lauter wird 
demgemaB der Ruf nach einer „dringend not- 
wendigen Regeometrisierung 66 (Schupp) des ge- 
samten Mathematikunterrichts. 


Mit Hilfe von DERIVE und der Unterstiitzung 
von Cabri-Geometre, GEOLOG oder EUKLID 1 
soil im folgenden ein Problem behandelt wer- 
den, das einerseits zur ,,Bekampfung 66 der For- 
menarmut im Geometrieunterricht geeignet ist, 
andererseits aber auch eine Briicke zur Algebra 
und Analytischen Geometrie bildet. Eine derar- 
tig „integrierende Behandlung 66 eines mathema- 
tischen Problems ist fur den Unterricht umso 
interessanter und notwendiger, als demnachst in 
Form des TI-92 ein Taschenrechner zur Verfu- 
gung steht, der u.a. Derive und Cabri-Geometre 
integriert hat. 

1. Introduction 

Traditional geometry teaching mainly deals 
with polygons and circles. Working with 
mappings is restricted on similarities which 
transform lines on lines and circles on cir- 
cles. Recently there are some attempts - 
influenced by software products like 
CABRI-Geometre - to involve at least the 
inversion on a circle. 

Beside the lack of shapes in teaching ge- 
ometry the separation in “Geometry” and 
“Algebra” in textbooks and exercise books 
will reinforce the pupils' impression of ge- 
ometry and algebra as two very separated 
mathematical fields. 

Using DERIVE - supported by CABRI, GE- 
OLOG or EUKLID 1 we will deal at this place 
with a problem which is able to bring new 
forms into Geometry teaching and to build 
up a bridge between Algebra and Geome- 
try as well. The new TI-92 makes DERIVE 
and Cabri available in one calculator and 
allows to work parallel on a geometric and 
on an analytic level as well. 

1 EUKLID is a shareware program offering the same 
features as CABRI and is running under WINDOWS. 
EUKLID can be downloaded from 
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2. Eine Verallgemeinerung der Punktspiege- 
lung 2 

Wie jeder Schuler in der Sekundarstufe lernt, ist 
die Punktspiegelung eine geraden-, strecken-, 
winkel- und kreistreue Abbildung. Durch eine 
kleine Modifikation der Abbildungsvorschrift 
lassen sich aber erstaunliche Veranderungen er- 
zielen. Nimmt man einmal Ernst, dass man 
niemals einen Punkt (im mathematisch ideali- 
sierten Sinn) zeichnen oder sehen kann, konnte 
man versuchen, einen Punkt durch eine Kreis- 
scheibe darzustellen und eine „Spiegelung“ an 
der Kreislinie durchzufiihren. 

Dieser Idee folgend, kommt man zu folgender 

Abbildungsvorschrift: 

Gegeben ist der Kreis k(M,r) und ein Punkt P. 

P wird am Schnittpunkt Z MP n k(M,r ) der 
Geraden MP mit der Kreislinie (punkt-) gespie- 
gelt und man erhalt den Bildpunkt P\ 


2. A Generalization of a reflection wrt a 
point 2 

As every secondary level student learns, 
the reflection wrt a point is a line-, seg- 
ment-, angle-, and circle true mapping. A 
little change of the mapping results in re- 
markable changes. We represent the point 
by a circle and then perform a “reflection” 
with respect to this circle. 

Following this idea we find the instruction 
for the respective mapping: 

Given is the circle k(M,r) and a point P. 

P is reflected at Z = MP n k(M,r) wrt to the 
circumference giving point P'. 

You can also reflect wrt to the tangent in Z. 

2 The mapping was discovered independently by me 
and Mr Wiesinger. We used it in teaching geometry. 
Mr Wiesinger will publish his experiences in the next 
future. 



(I used Tl-N Spire to reproduce Thomas' original CABRI-screen shots. You can do it easily on 
the Tl-handhelds TI-92, Voyage 200, too. You can also demonstrate the mapping with DE- 
RIVE using slider bars. Josef) 
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2 2 2 
#1: (x - m) + (y - n) = r 

#2: [0 := [m, n], P := [u, v]] 

v - n 

#3: y - v = ■ (x - u) 

u - m 

u - m 

#4: y - (x - (S(P)) ) 

2 v - n 1 

Sliders for m,n,r, and u,v. 

P* 

#5: S(pt) := 0 + 

|pt 

#6: P_(pt) := 2-SCpt) 

#7: S(P) 

#S: P_(P) 

3. Abbildungseigenschaften / Properties of this Mapping 

Um die Abbildung kennenzulernen, wird man zunachst Phanomene untersuchen, die sich beim Abbil- 
den von Punkten ergeben. Im vorliegenden Fall beobachtet man z.B.: 

• Die Kreislinie k(M,r) ist ein Fixpunktkreis. 

The circle k(M,r) is a fix circle. 

• Die Kreislinie k(M,2r) wird auf den Kreismittelpunkt M abgebildet. 

The image of the circle k(M,2r) is the center M. 

• Die Kreislinie k(M,Ar) wird auf die Kreislinie k(M,2r) abgebildet. 

The image of the circle k(M,4r) is the circle k(M,2r). 

• Niemals liegen P und P r gleichzeitig innerhalb des Kreises k(M,r). 

P and P ' can never lie within the circumference of k(M,r) at the same time. 

r :e Real (0, ®) 

P_C[m + r-COSCoc), n + r-SIN(oc)]) = [r-COS(oc) + m, r-SIN(oc) + n] 

P_([m + 2'r'COS(oc), n + 2 ■ r ■ SIN(oc) ]) = [m, n] 

P_([m + 4'r'COS(oc), n + 4 ■ r ■ SIN(oc) ]) = [m - Z-r-COSCoc), n - 2-r-SIN(oc)] 

Neben diesen grundlegenden Phanomenen, die sich leicht erklaren lassen, liefert die Untersuchung des 
„Symmetrieverhaltens“ interessantere Beobachtungen. Zur Untersuchung spezieller Symmetrien wie 
Dreh-, Verschiebungs-, Achsen- oder Punktsymmetrie eignet sich die Betrachtung des Bildes, das sich 
beim Abbilden einer Standardfigur, etwa eines groBen „F“ ergibt. Die oben genannten Geome- 
trieprogramme erlauben es, ein Konstruktionsobjekt (etwa den Urpunkt P) zu bewegen und gleichzei- 
tig die „Ortslinie“ eines konstruktiv abhangigen Punktes (etwa den Bildpunkt P *) zu protokollieren. 

We want to investigate special symmetries and observe the mapping of a standard figure - 
eg an uppercase “F“. The geometry programs mentioned above allow to move the original 
point P and trace the locus of its picture point P\ 


- o 

■ r 

- o| 

- pt 



In den obigen Abbildungen wurde P entlang 
des groBen „F“ bewegt; die Punktwolke ist die 
protokollierte Ortslinie des Bildpunkts P' . Das 
verzerrte Bild der linken Abbildung lasst keine 
besondere Symmetrieeigenschaft erkennen. 
Dies andert sich, wenn man sich mit der Urfigur 
von der Kreislinie entfemt oder, was gleichbe- 
deutend ist, den Kreisradius sehr klein macht. 
Plotzlich verhalt sich die Abbildung ganz an- 
ders und liefert ein Bild, das man von der 
„normalen“ Punktspiegelung am Punkt M er- 
warten wiirde. 

Bei der vorliegenden Abbildung handelt es sich 
also in der Tat um eine Modifikation einer 
Punktspiegelung, bei der im Grenzfall fur groBe 
Entfemungen noch die urspriinglichen Eigen- 
schaften erkennbar werden. 


In both pictures above point P was moved 
along the “F”. We can see the locus of the 
mappings P\ One cannot recognize any 
symmetry in the distorted left picture. But 
things change by moving the base figure 
away from the circle, or equivalently de- 
creasing the radius of the circle. Suddenly 
the mapping produces a picture which 
could be expected as the result of an “ordi- 
nary” reflection wrt to point M. 

So this mapping is indeed a modification of 
a reflection wrt to a point which for large 
distances shows the original properties. 


This is the DERIVE realisation using the sliders simulating a dynamic geometry program - 
but based on the mathematical models of the original projects and the mapping (see expres- 
sions #5 and #6 from above). 

The TABLE command allows to plot the loci for the three parts of the “F”. Moving the sliders 
for m and n changes the position of the center of the circle while the slider for r varies the ra- 
dius of the circle. 


F := [[-4, -2] + t ■ [0, 6], [-4, 4] + t-[3, 0], [-4, 2] + t ■ [2 , 0]] 
j-p-CF^), p_(F_y p -oy] 

[ (TABLE (P_(F ), t, 0, 1, 0.02))u2, (TABLE(P_(F ), t, 0, 1, 0.02))u2, CTABLE(P_(F ), t, 0, 1, 0.02))u2 
L 1 2 3 
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Diese ersten Beobachtungen konnen in der Sekundarstufe I verschiedenen Zwecken dienen: 

Zum einen konnte man die Abbildung vor der Achsen- und Punktspiegelung behandeln, um letztere 
dann als Grenzfalle (r 0) zu gewinnen. Bei diesem Vorgehen kann die „Emfachheit“ der Abbil- 
dungseigenschaften der Kongruenzabbildungen als Erleichterung empfunden werden - und nicht wie 
bei der herkommlichen Behandlung als „selbstverstandlich“. 

Zum anderen konnte die Abbildung am Ende einer Unterrichtssequenz zu Ahnlichkeitsabbildungen 
stehen, um wenigstens einmal eine nichtgeradentreue Abbildung vorzustellen oder vom Schuler pha- 
nomenologisch untersuchen zu lassen. 

One could use this mapping at the end of a teaching unit about similarity mappings in order 
to present the students at least one non line-true mapping and let them investigate this map- 
ping. 


4. Die Bilder von Geraden / The images of lines 

Nach dem Verhalten von Punkten unter der Abbildung wird man sich im nachsten Schritt um die Bil- 
der von Geraden bemuhen. Dazu variiert man zunachst einen Urpunkt auf einer Geraden und betrach- 
tet die entstehende Ortslinie des Bildpunktes P\ Fur verschiedene Abstande der Urgeraden vom 
Kreismittelpunkt sind in den folgenden Abbildungen einige Geradenbilder dargestellt. 



wwzwTm i 

- ^ <[\ 

t — - 

1 

MAIN I 1 EG All’ll 

u — — 
p 

FUNC 







m ' FUNC 

H^h ^ DDI i 




(Compare the TI-92 / Voyage 200 screens with the graphs produced with CABRI or EUKLID. 
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For different distances of the line containing the moving point P one obtains different curves 
for the trace of P (loci). We can find out that there is an asymptote. We try to find an alge- 
braic representation of these curves in order to investigate them more accurately. 


Wie die Bildfiguren erkennen lassen, kann man zunachst eine grobe Einteilung in Kurven mit „Schlei- 
fen“, einer „Spitze“ und mit „Beulen“ treffen. Genauere Untersuchungen deuten darauf hin, dass die 
Bildkurven eine Gerade als asymptotische Naherungskurve besitzen, was sich anschaulich auch be- 
griinden lasst, wenn man die oben gemachten Beobachtungen beziiglich des „Symmetrieverhaltens“ 
beachtet. 

Von den dargestellten Kurven mochte man nun eine mathematische, d.h. algebraische Darstellung er- 
halten, um sie naher untersuchen zu konnen. Dazu definieren wir ein Koordinatensystem mit dem 
Kreismittelpunkt als Ursprung (dieser bietet sich als einziger ausgezeichneter Punkt der Abbildung an) 
und legen die Achsen „kanonisch“, d.h. die x-Achse von links nach rechts und die y-Achse von unten 
nach oben. 


As we can recognize from the image curves one can have a rough classification: curves with 
a “loop”, curves with a “vertex” and curves with a “bump”. We have the idea that there is a 
straight line as asymptote (notice the symmetry behaviour). 

We want to find an analytical presentation of the curves. For this purpose we embed the fig- 
ure in a system of coodinates with the centre of the circle as origin and the axes as usual. 


Geben wir dem Kreis den Radius 1, so hat ein Kreis- 


^coscr^ 


. Ein 


ysma j 


punkt Z die Koordinatendarstellung z = 

Punkt P mit der Entfernung r vom Ursprung hat die 

Fur den Bildpunkt P" er- 

^(2-r)cos6r^ 


f rcosa ^ 


v rsm a j 


Darstellung p = 
halt man nun p = p + 2 (z - p) = 


(2 - r) sin a 


x 

Ersetzt man in dieser Darstellung cos a = — , 

r 

sin a = — und r = yj x 2 + y 2 , dann erhalt man fur 
r 

den Bildpunkt 


— r 

p = 


x^2 -yjx 2 + y 2 j 

^( 2 - 7777 ) 

v V^ 2 + ~y 2 y 


#1: 

k := [COS (a) , SIN (a) ] 


#2: 

P := [ 

r-CQSCoO, r-SIN(a)] 


#3 : 

pi := 1 

p + 2- (-p + k) 


#4: 

pi : = 

[(2 - r) -COSta), (2 - r)-SIN(a)] 

#5 : 

x = r 

■ COS(a) 


#6: 

y = r 

■SIN (a) 


#7: 

Pi = 

[(2 - r) -COSta), (2 - r)-SINta)] 

#8 : 

Pi = 

' (2 - r) -x (2 - r)-y ' 



r r 




You can follow the calculation done with and 
without DERIVE for obtaining the coordinates of 
the image of point P = P\ 
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Um das Bild einer Parallelen zur y-Achse zu erhalten substituiert man in dem oben erhaltenen mar- 
kierten Ausdruck x durch einen festen Wert t (= Abstand von der y-Achse) und y durch einen Kurven- 
parameter a. Im Bild auf der vorigen Seite sind die Kurven fiir ^-Werte von -3 bis -0,5 mit der Schritt- 
weite 0,5 dargestellt (mit Hilfe des VECTOR-Befehls). 

For obtaining the image of a vertical line one has to substitute in the above highlighted ex- 
pression x by a constant value t (= distance from the y-axis) and y by the curve parameter a. 
The graph shows the family of curves for -3 < t < -0.5 with an increment of 0.5 (using the 
VECTOR-command). 

Eliminiert man aus der Parameterdarstellung eines Geradenbildes 

t[l-y]t 2 +a 2 j a{l-\lt 2 +a 2 j 

\lt 2 +a 2 yjt 2 +a 2 

V ) 

den Kurvenparameter a, so erhalt man als Beziehung zwischen den Koordinaten x und y der Bildkur- 
ven die algebraische Gleichung (x 2 + y 2 )(x + ^) 2 =4x 2 , also die Konchoidengleichung (vgl. Folge 5 

des Kurvenlexikons in DNL#16). In der folgenen Abbildung sind die einzelnen Schritte mit DERIVE 
durchgefuhrt. 


We eliminate parameter a obtaining the algebraic equation (x 2 + y 2 )(>x + /) 2 =4x 2 for the 
family of image curves. These are Conchoids (see Lexicon of Curves (5) in DNL#16). 


#12: 


2 2 2 2 
t'(2 - A t + a )) a ■ ( 2 - At + a )) 

* = 1 y = 

2 2 
At + a ) 


#13: 


2 2 


At + a ) 


2 

2 

t ■ C 2 - At + 

« )) 

2 2 


At + a ) 


2 

2 

ot'(2 - + 

a )) 


2 2 
At + a ) 


#14: 


x t 


y ct 


#15 : GUEST 


2 2 

t ■ C 2 - VCt + a )) 

t-y 


f 2-SIGNCt)' \x I ^ 

x ~ t 

2 2 
At + a ) 

X 


2 2 

^ A* + y ) 


#16: 


x + t 


2-SIGNft) ■ | a | 

2 2 
A* + y ) 


2 4-x 

#17: O + t) = 

2 2 
x + y 

2 2 2 2 

(_x + t) ■ (a + y ) = 4 'X 


#18: 
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5. Abschliefiende Bemerkungen / Final Comments 


Mit der Erkenntnis, dass es sich bei den betrachteten Geradenbildern um Konchoiden des Nikomedes 
handelt sind fur die Sekundarstufe II neue mathematische „Turen“ aufgestoBen. Es stellen sich etwa 
die Fragen: 


• Wozu braucht man eine Konchoide (vgl. Folge 5)? 

• Gibt es andere Konchoidenkonstruktionen? 

• Wie sieht die urspriingliche Konchoidenkonstruktion aus? 

• Welche Kurven ergeben sich, wenn man Kreise abbildet? 

• Wie sehen die Bilder von Ellipsen, Parabeln, Hyperbeln aus? 

• Wie lassen sich die Unterschiede zwischen den „Parameterplots“ und den „Implizit-Plots“ 
bei DERIVE erklaren (vgl. unten stehende Abbildung)? 

• Kann man die Konstruktionsvorschrift so erweitern, dass die Bilder mit denen der implizi- 
ten Plots ubereinstimmen? 




Die Auswahl der Fragen soil nur einen kleinen Hinweis geben, wie das geometrische Ausgangsprob- 
lem sowohl algebraisch, analytisch und elemetargeometrisch weiter durchdrungen werden kann. Die 
Beantwortung dieser Fragen wiirde den Rahmen dieses Artikels sprengen und sei dem interessierten 
Leser iiberlassen. Fur Antworten ware ich sehr dankbar. 


Learning that the images of straight lines are Conchoids new “mathematical doors” for sec- 
ondary level 2 are pushed open. A couple of new questions can be posed, as 

• For which purpose do we need a Conchoid (see Lexicon 5)? 

• Are there other consructions for this curve? 

• How does the original construction of this curve look like? 

• What are the images of circles, ellipses, hyperbolas, parabolas, ...? 

• How can one explain the difference between the parameter plot and the implicit 
plot in DERIVE (see figures above)? 

• Is it possible to extend the construction in such a way that both plots coincide? 

I really would appreciate any answers of members of the DERIVE community 


6. Literatur / References 

Bock/Werge, Finden von Vermutungen durch funktionale Betrachtungen, MidS 31 (1993) 2 
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In 1996 I didn't have time and possibly not the means to react on Thomas Weth's invitation 
and challenges. One of my favourite DERIVE features - the slider bars - inspired me to im- 
mediately try mapping circles and hyperbolas. 

The procedure is very easy: one has only to substitute x and y in expression #9 by the re- 
spective parameter representation of the curves to be mapped. 

I didn't try to find implicit forms of the resulting curves - and I doubt if this is possible. 

Josef 


Images of circles [(m,0);r]: 


2 2 2 2 
*.(2 -VO + y )) y-c 2 - VC* + y )) 


2 2 2 2 
VO + y ) VO + y ) 


#20: [m + r -COS (a) , r -5IN(cO] 


2 2 2 2 
(m + r 'C05 (h)) ■ (_2 - VCO + r .C05(a)) + ■ SINCa)) )) O ■ SINCa)) ■ (2 - VC Cm + r .CQS(a)) + (r.SIN(a)) )) 


2 2 
VC Cm + r-CGSCcO) + Cr-SINCa)) ) 


2 2 
VC Cm + r-COSCct)) + Cr-SINCa)) ) 


2 2 2 2 
Cr-C05Ca) + m)-0 - VO -m ■ r -C05W + m + r )) r ■ SINCa) -O - VO -m ■ r -COSCa) + m + r )) 


2 2 

VO 'in ■ r -CGSO) + m + r ) 


2 2 

VO 'in ■ r -COSO) + m + r ) 



Images (red) of one branch of a hyperbola (black): 

#23: [m + a-CQSHCa), b-SINHCa)] 


#24: 


2 2 

Cm + a-CGSHCa)M2 - VCO + a-COSHCa)) + Cb-SINHCa)) )) 


2 2 
VCC«i + a 'COSHCa)) + C^-SINHCa)) ) 


2 2 

Cb-SINHCa )) -(2 - VC Cm + a'CQSHCa)) + Cb-SINHCa)) )) 


2 2 
VC Cm + a.CQSHCa)) + Cb-SINHCa)) ) 
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JULIA Sets with DERIVE 

Hartmut Kiimmel, Biedenkopf, Germany 

This is a short - self explanatory - file for producing and plotting Julia sets. Hartmut has also submit- 
ted a paper “The representation of plane figures with DERIVE”. Josef 

— COMPLEX. dfw - A Tool for Complex Numbers — 

— Tools for the Escape-Algorithm — 

#1: [maxi ter := 32, maxdist :=, c :=] 

2 

#2: Q(z, c) := z + c 

ITNUMBER(z , it) := 

If it > maxi ter v ABS(z) > maxdist 
#3: it 

ITNUMBER(Q(z , c) , 1 + it) 

#4: ESC_NUMB (start) := ITNUMBER(start , 1) 

#5: ESCAPE(start , c) := ITERATES (Q(z , c), z, start, 10) 

POINTS(v) := VECTOR([RE(v ), IM(v )], j, DIMENSION(v)) 

#6: L j j J 

First Example: 

#7: [maxdist := 2.5, c := -0.567 + 0.456-i., a := 0.325 + 0.325 -l] 

#8: ESC_NUMB(a) = 10 

#9: ESCAPE(a, c) 

#10: POINTS (ESCAPE (a, c)) 



— Tools to represent Julia Sets — 

#11: GRID(xst, xend, dx, yst, yend, dy) := APPEND(VECTOR(GRIDL(xst, xend, dx, 

y), y, yst, yend, dy)) 

#12: CATCH(pts) := SELECT(ESC_NUMB(z) > maxi ter, z, pts) 

#13: GRIDL(xst, xend, dx, y) := VECT0R(x + y -i, x, xst, xend, dx) 
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Example: 

#14: CATCH ([0.325 - 0.325-l, 0.5 - 0.5-l, 0.75 - 0.75-l]) 

#15: [0.325 - 0.325-l, 0.5 - 0.5-l] 

#16: testl := GRID(-2, 2, 0.1, -1.5, 1.5, 0.1) 

#17: POINTS(testl) 

Plot #17: points Size Small and Color Grey 

#18: POINTS (CATCH (testl)) 


Plot #18: points Size Medium and Color Red 



Two more examples: 



#19: 

c := -1 

#22: 

c := 0.32 + 0.043 - l 

#20: 

POINTS(testl) 

#23: 

POINTS(testl) 

#21: 

POINTS (CATCH (testl)) 

#24: 

POINTS (CATCH (testl)) 
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Jan Vermeylen from Kapellen, Belgium, found another approach to represent Julia sets: 


#1: 

Notation := Decimal 


#2: 

NotationDigi ts := 3 


#3: 

resol := 750 



#4: 

IT(c) := ITERATES((2 • RAND0M(2) - l)-/(z - 

c) , z, 0, resol) 

#5: 

BEELD(z) := [RE(z), IM(z)] 


#6: 

FIG(p) := VECTOR(BEELD(ELEMENT(p, n)), n, 

10, resol) 

#7: 

DULIA(c) := FIG(ITCc)) 


#8: 

( 3 1 
DULIA 
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#11: DULIAC0.1 - 0.8-L) 



“Fatou Dust” 
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Investigate Distributions of 
Functions of Random Variables 

Josef Bohm, Wiirmla, Austria 

In 1996 DERIVE had no real programming features implemented. We had to work with 
functions calling other functions which was sometimes very complicated. We had to use the 
ELEMENT-function instead of SUB = i and there were some other uncomfortable facts 
compared with DERIVE 5 or DERIVE 6. So I decided to rewrite the DERIVE code and the 
respective possible handout for the students. Josef 

This is the handout for my students for trying finding the relation between expected values and vari- 
ances of discrete random variables X, 7, Z, ... and the expected values and variance of random vari- 
abless U, V, W, ... which are functions (preferable linear) of X, 7, Z, . . . 

Investigate Distributions of 
Functions of Random Variables 

Assume we know the distributions of three random variables X, Y and Z: 


X 

0 2 3 4 

Y 

3 4 5 

Z 

-1 1 

p(X) 

0.1 0.2 0.3 0.4 

P(Y) 

0.2 0.3 0.5 

P(Z) 

0.8 0.2 


and we would like to find the distribution of random variables U. !'. IV, ... which are functions of X, Y 
and Z. These distribution tables can be used to calculate the means and variances of the new variables. 
Our aim is to find out if there are some relations between the means and variances of the random vari- 
ables and the random function values. 

Load the file LINRAND.MTH as an Utility file. 

Hint: Edit the distribution tables as matrices, eg: vl := [0,0.1;2,0.2;3,0.3;4,0.4], 

Then investigate the distributions of: 

U=3X 

V = 3X- 2 
W= 3X+2Y-4 

T=2X+3Y-4Z-2 

S= -4X- 5Y 

Q = 2X 2 

by calculating their expected values and variances. 

Try to find the relation between expected values and variances of the given random variables and the 
newly created random variables. If you cannot find a satisfying answer or if you would like to confirm 
your conjectures then generalize, e.g. U= aX, V= aX+ fi, W = aX-fi 7-y, etc. 
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pr_fc(dists, function, variables) returns the "condensed" table, i.e. without showing the 
values for the single variables: 

' -5 0.08 

-13 0.16 

-17 0.24 

-21 0.32 

#5: pr_fc([vl, v3], M, [x, z]) = 

15 0.02 

7 0.04 

3 0.06 

. -1 0.08 


You may prefer sorting the table: 
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' -21 0 . 32 
-17 0.24 

-13 0.16 

-5 0.08 

#6: SORT(pr_fc([vl, v3], M, [x, z])) = 

-1 0.08 
3 0.06 

7 0.04 

.15 0.02 

If there are more than 2 variables it can happen that some function values appear 
more often than only once. Then pr_ff(condensed and sorted distribution table) 
will return a new table with collected (= added) probabilities of same function values. 

See one example: 

' -l 0.512 
1 0.128 
-3 0.128 

-1 0.032 

#7: pr_fc([v3, v3, v3], x - y + z, [x, y, z]) = 

1 0.128 
3 0.032 

-1 0.032 

. 1 0.008 

' -3 
-1 

#8: pr_ff (S0RT(pr_fc([v3 , v3, v3], x - y + z, [x, y, z]))) = 

1 

. 3 

We proceed with variable M and calculate its mean and variance. 
exp_val (condensed distribution table) gives the mean and 
vari (condensed distribution table) gives the variance: 



exp_val (pr_fc([vl, v3], M, [x, z])) 
#9: 

. vari (pr_fc([vl, v3], M, [x, z])) 

Mean and variance of vl and v3 are: 



- 12.6 ' 
87.84 . 


# 11 : 


0.64 
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I repeat from above: function M(X,Z) = -4X + 10Z +5. 

Do you find a relationship between the means and variances? 

Start with one variable functions then it will be easier!! Generalize!! 


I printed the file in 1996, so I am doing in 2010, too. If you have DNL#22 available then you 
are invited to compare! 

pr_f(dists, f, vs, dl, d2, dB, dt) := 

Prog 

If DIM(dists) = 1 
Prog 

dl := di sts|l 

dt VECTOR([SUBST(f , vs|l, dl|i|l), dl|u2], i, DIM(dl)) 

If DIM(dists) = 2 
Prog 

dl := di sts|l 
#1: d2 := dists|2 

dt APPEND(VECTOR(VECTOR([dl|ia, d2ija, SUBST(SUBST(f , vs|l, 
dl|i|l), vsi2, d2|j|l), dl|i i2*d2|j|2] , i, DIM(dl)), j, 

DIM(d2))) 

If DIM(dists) = 3 
Prog 

dl := dists^l 
d2 := dists|2 
d3 := dists|3 

dt APPEND(APPEND(VECTOR(VECTOR(VECTOR([dl|i|l, d2|j|l, d3|k|l, 
SUBST(SUBST(SUBST(f , vs|l, dl|i|l), vs|2, d2ija), vs|3, 
d3|k|l), dl|i |2-d2|j|2-d3|k|2] , i, DIM(dl)), j, DIM(d2)), k, 
DIM(d3)))) 
dt 

pr_fc(dists, f, vs, dl, d2, d3, dt) 

If DIM(dists) = 1 

#2: pr_f(dists, f, vs, dl, d2, d3, dt) 

(pr_f(dists, f, vs, dl, d2, d3, dt))u[DIM(vs) + 1, DIM(vs) + 2] 

pr_ff(d, f, dt) 

Prog 

dt := [d^l] 

Loop 

d := REST(d) 

If d = [] 

#3: RETURN dt 

f := FIRST(d) 

If fil - (FIRST(REVERSE(dt)))|l 
Prog 

dt := APPEND(REVERSE(REST(REVERSE(dt))) , 

[ [ ( FI RST( REVERSE (dt))) 4,1, (FIRST(REVERSE(dt ))) \2 + f|2]]) 

Prog 

dt := APPEND(dt, [f ]) 


#4: 

exp_val (dt) 

dtj,j,l*dt,u2 

#5: 

vari (dt) := 

fDIM(dt) 2 T 2 

I I dt -dt - exp_va1 (dt) 

l i=l i , 1 i , 2 J 

#6: 

Notation := 

Decimal 
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Finding a Limit via 
Geometric Reasoning 

Carl Leinbach and Marvin Brubaker, USA 

Before we begin this investigation, adjust the graphics window to our needs using the Set > As- 
pect Ratio >1:1 option. The screen should look similar to the figure below. 


m 


(Zl[§j[X| ES 2D plot 1:1 




• 0.2 


■ 0.2 


0.2 0.4 0.6 0.8 1 


- 0.2 


Consider the following sequence of points: P n = < 


[0,0] 

n- 0 

[0,1] 

n - 1 

[1,0] 

n-2 

\Pn - 3 +lPn - 2 

otherwise 


Notice that this sequence is defined recursively. DERIVE allows us to make recursive definitions. We 
use the IF statement. 


P (n) : = IF (n=0, [0, 0] , IF (n=l, [0,1] , IF (n=2, [1, 0] , 1/2P (n-3) +1/2P (n-2) ) ) ) 

In this case we had to nest the IF statements three deep. That is because we had three special cases. 

This function, because of its recursive nature, is slow to evaluate for an n of any size, whatsoever. 
Nonetheless, author 

VECTOR ( P (n) ,n, 0, 10) 

and plot the sequence. 

(It is not necessary to simplify the expression - giving a matrix of points. But take care that you have 
activated the Option > Simplify before plotting or Approximate before plotting in the plot window. 

Set the Points in the Display Options Connected and Size Small. 

The next figures show the evaluation of the first 10 terms of the sequence and also the first 20 terms. If 
we move the crosshair on the graph where the plot is dense, i.e., the point of apparent convergence we 
get a reading of approximately [0.4, 0.4]. 
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We can zoom in and the we read off the coordi- 
nates of the crosshair [0.40029, 0.40042], 

We can show the last term of the sequence given 
right above and we get a similar result: 
[0.40039..., 0.40039 ...]. 

Of course, we had not proved any result. How- 
ever, the visual evidence is convincing that a 
limit does exist ([0.4, 0.4]?) and we have a visual 
illustration of the process of convergence. 



FIRST(REVERSE (VECTOR (P(n) , 
FIRST(REVERSE (VECTOR (P(n) , 


0 , 20 ))) 


205 205 


512 512 

0 , 20 ))) = [ 0.400390625 , 0 . 400390625 ] 


As Carl wrote, the recursive function is slow - try for n = 50! With DERIVE 5 and higher we 
can write a small program - without applying the interesting recursive function from above - 
which allows to calculate much more elements of this sequence. 


pts(n, pt) 

Pnog 

pt :z [0, 0; 0, 1; 1, 0] 

k::4 
#14: Loop 

If k > n 
RETURN pt 

pt := APPEND(pt, [1/2- (pt|(k - 3) + pt|(k - 2))]) 
k :+ 1 


#1B : 

pts(50) 




#16: 

FIRST(REVERSE(pts(50))) = 

" 33B5443 

6710887 " 


. S3SS60S 1 16777216 . 


#17: 

FIRST(REVERSE(pts(50))) = 

[0.3999999761 

, 0,4000000357] 


The challenge is still there: Proof that the limit is [0.4, 0.4]! 
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DREIECK.MTH - TRIANGLE.MTH (1) 

Berhard Wadsack, Vienna, Austria 


Maybe that you are remembering Berhard 's nice report about his lesson using a 10m DE- 
RIVE print out of n in a Viennese grammar school in DNL#17. Some time ago Bernhard 
submitted a paper to calculate and to plot the “remarkable points” (die merkwurdigen Punkte) 
of a triangle: orthocenter (Hohenschnittpunkt), circumcenter (Umkreismittelpunkt), incenter 
(Inkreismittelpunkt) and centroid (Schwerpunkt). In the following DERIVE file Berhard refers 
to David Sjostrand's contribution “CAS and Spreadsheets” from DNL#13. He does not derive 
the coordinates of circumcenter and incenter but he uses David's results. Berhard wrote that 
this contribution gave the impetus for his work which kept him busy for months. 

I like this contribution because it can be a possibility to encourage modular working. Groups 
of pupils should investigate the properties of a triangle and use their knowledge to produce a 
“Black Box” for further use. The file provides the derivation of nearly all formulae - excep- 
tions see above - and the results are collected in several numerical tables and in some ex- 
pressions which can be plotted immediately. 

I don't change Bernhard's original variables' names, but I try to give a translation and in- 
clude the English forms in the file. So you can load TRIANGLE.MTH as a utility file and start 
working. You can find worked examples on page 37. 


The expressions giving numerical results in tables : 

bp op Beriihrpunkte des Inkreises / Osculation pts of the incircle 

fusspunkte pedpoints HohenfuBpunkte / Pedal points of the altitudes 

The expressions for immediate plotting : 


seiten 

sides 

Dreieck / Triangle 

0 < p < 1 

umkreismpkt 

circcenter 

Umkreismittelpunkt / Circumcenter 


umkreiscirccirc 

Umkreis / Circumcircle 0 < p < 2n 

seitsymm 

perpbisecs 

Seitensymmetralen / Perp. bisectors 

-5 < p < 5 

inkreismpkt 

incenter 

Inkreismittelpunkt / Incenter 


inkreis 

incircle 

Inkreis / Incircle 

0 < p < 2n 

winkelsymm 

angbiss 

Winkelsymmetralen / Angle bisectors 

0 < p < 1 

hoehenschnpkt 

orthocenter 



hoehen 

altitudes 

Hohen / Altitudes 

-5 < p < 5 

schwerpunkt 

centroid 



schwerlinien 

medians 

Schwerlinien / Medians 

0 < p < 1 


The list will be accomplished in DNL#23 offering TRIANGLE.MTH (2) 

Comment from 2010: I didn't change very much and left the file in its original form as far as possible. 
I included the English variable names. The file can be used by English and German speaking users as 
well now. At least I am hoping so, Josef 
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#2B : 
#26: 
#27: 
#28: 
#29: 
#30: 
#31: 
#32: 
#33: 


#36: 

#37: 

#38: 

#39: 

#40: 

#41: 

#42: 


#44: 
#4B : 
#46: 
#47: 
#48: 
#49: 
#S0 : 

#S1: 


InputMode := Word 

[xl :=, yl :=, x2 :=, y2 :=, x3 :=, y3 :=] 

[pi := [xl, yl], p2 := [x2 , y2], p3 := [x3 , y3]] 

[eckpunkte := [pi, p2 , p3], vertices := [pi, p2 , p3]] 

Die Parameterdarstellungen der 3 Sen ten / parameter form of the sides 
:= [x2 + t-(x3 - x2) , y2 + t- (y3 - y2)] 

b_ := [x3 + t- (xl - x3) , y3 + t- (yl - y3)] 

c_ := [xl + t- (x2 - xl) , yl + t- (y2 - yl) ] 

[sei'ten := [a_ , b_, c_], sides := [a_, b_ , c_]] 

Andere Hogli'chkeit das Drei'eck zu zei'chnen / other way to plot the triangle 
DREIECK(xl, yl, x2, y2, x3, y3) 


TRIANGLE (xl, yl, x2, y2, x3, y3) 

Zum Umkrei'snn' ttelpunkt / Center of the ci' rcumci rcl e : 


Umkrei's radius / its radius: 

2 2 

r := [ (xu - xl) + (yu - yl) ) 

Berechnung der Koondinaten von U(xu,yu): 

22 2222222 
xl -x2 - xl -x3 - xl-x2 - xl-y2 + xl-x3 + xl-y3 + yl -x2 - yl -x3 + x2 -x3 - x! 

yu := 

2- (xl-y2 - xl-y3 - yl-x2 + yl-x3 + x2-y3 - y2-x3) 

2222 222222 
xl ■ y2 - xl -y3 + yl -y2 - yl -y3 - yl-x2 - yl-y2 + yl-x3 + yl-y3 + x2 -y3 + y2 ■ 

xu := 

2-(xl-y2 - xl-y3 - yl-x2 + yl-x3 + x2-y3 - y2-x3) 

u_ := [xu , yu] 

[umkrei'smpkt := u_, ci rcci rccenter := u_] 

Parametendarstellung des Umkrei'ses / parameter form of the ci rcumci rcl e 
[umkrei's := [xu + r-COS(t), yu + r-SIN(t)], ci rccenter := umkrei's] 

Die Streckensymmetralen halbi'eren die Sei'ten des Drei'ecks und stehen auf di'ese normal 
Di e Sei' tenmi' tten si'nd m_a, m_b und m_c . 

The perpendicular bisectors of the sides pass the midpoints m_a, m_b and m_c of the sides. 

p2 + p3 pi + p3 pi + p2 

m_a := , m_b := , m_c := 

2 2 2 
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r u 

f u 1 


r u t 

( u 

— . 

la 


lb ■ 

lc 

L 2 

L 2 J 


L 2 J 

L 2 JJ 


#84: Inkrei's radius / radius of inci'rcle: 

2-f 


#86: Inkrei's / inci'rcle: 

#87: [inkrei's := [x_i + p-CQS(t), y_i' + p-SIN(t)], inci'rcle := inkrei's] 

#88: Die Beruhrpunkte des Inkrei'ses / the osculation points of the inci'rcle ta, tb, tc : 

#89: [nl := [x_i + t- (y3 - y2) , y_i' + t- (x2 - x3)], n2 := [x_i + t- (yl - y3) , y_i + t- (x3 - xl)]] 

#90: n3 := [x_i + t- (y2 - yl) , y_i' + t- (xl - x2)] 

p. (y3 - y2) P- (x2 - x3) 

#91: ta := x_i + , y_i + 

la la 


p. (yl - y3) p- (x3 - xl) 

#92 : tb := x_i' + , y_i + 

lb lb 

r P- (y2 - yl) P- (xl - x2) 

#93 : tc := x_i' + , y_i + 

lc 1c 

#94: [ta, tb, tc] 

#9B: Beruhrpunkte des Inkrei'ses bp / osculation points of the inci'rcle op: 


tc := x_i 


I" Beruhrpunkte des Inkrei'ses [x, y] 1 I" osculation pts of inci'rcle [x , y] T 

an der Sei'te a: ta on side a a: ta 

#96 : bp := , op : = 

an der Sei'te b: tb on side b b: tb 

L an der Sei'te c: tc J [ on side c: tc J. 

#97 : 6464464646446464464646446464464646446464464646446464464 

#98: Zum Hohenschm' ttpunkt / the Orthocenter: 

#99: Parametendarstellung der 3 Hohen / parameter form of the altitudes: 

#100: [h_b := [x2 + t2-(yl - y3) , y2 + t2-(x3 - xl) ] , h_a := [xl + tl- (y 3 - y2) , yl + tl- (x2 - x3)]] 
#101: h_c := [x3 + t3- (yl - y2) , y3 + t3-(x2 - xl) ] 

#102: [hoehen := [h_a, h_b, h_c], altitudes := hoehen] 

#103: Urn die Koondi'naten von H in ei'nem bel . Drei'eck zu berechnen, schnei'den wi' r h_a mi't h_b: 

#104: Wi'e intersect h_a and h_b to obtain the orthocenter H: 

#10S : h_a = h_b 

#106: S0LVE(tl-Cy3 - y2) + xl = t2-(yl - y3) t x2 a tl- (x2 - x3) + yl = t2- (x3 - xl) + y2, [tl, t2]) 


xl - xl- (x2 + x3) + x2-x3 + (yl - y2) ■ (yl - y3) 

tl = a t2 

xl-(y2 - y3) + x2- (y3 - yl) + x3-(yl - y2) 


xl- (x2 - x3) - x2 + x2-x3 + 
xl- Cy2 - y3) + x2-(y3 - y! 


SUB ST [tl- (y3 - y2) + xl, tl- (x2 - x3) + yl] , tl, 


xl- (y2 - y3) + x2-(y3 - yl) + x3- (yl - y! 






D 


-N-L#22 


B. Wadsack: 


DREIECK.MTH - TRIANGLE.MTH 


( 1 ) 


p35 


#109: 


2 2 2 

(x2 - 2-x2-x3 + x3 + (y2 - y3) ) ■ (yl - y2)-(yl - y3) x2- (yl - y2) + x3-(y3 - yl) 


Cy3 - y2)-(xl-(y2 - y3) + x2- (y3 - yl) + x3-(yl - y2)) y3 - y2 


3 2 2 2 2 2 

(x2 - 3-x2 -x3 + x2-(3-x3 + (y2 - y3) ) - x3- (x3 + (y2 - y3) ))-(yl - y2)-(yl - y3) 


2 

Cy2 - y3) ■ (xl- (y2 - y3) + x2-(y3 - yl) + x3- (yl - y2)) 

2 2 2 
x2 . (yl - y2) - x2-x3-(2-yl - y2 - y3) + x3 ■ (yl - y3) + yl- (y2 - y3) 


2 

(y2 - y3) 

#110: Factor #109 


# 111 : 


xl- (x2- (yl - y2) + x3- (y3 - yl)) + (y2 - y3)-(x2-x3 + (yl - y2) - (yl - y3)) 
xl-(y2 - y3) + x2- (y3 - yl) + x3-(yl - y2) 


2 2 2 2 2 
xl . (x2 - x3) - xl-(x2 - x3 + yl- (y3 - y2)) + x2 -x3 - x2-(x3 + y2- (yl - y3)) + x3-y: 


xl-(y2 - y3) + x2-(y3 - yl) + x3- (yl - y2) 


2 

xl - xl- (x2 + x3) + x2-x3 + (yl - y2)-(yl - y3) 

#112: xh := (y3 - y2) + xl 

xl-(y2 - y3) + x2-(y3 - yl) + x3- (yl - y2) 


#114: [hoehenschnpkt [xh, yh], orthocenter [xh, yh]] 

#11S : Fusspunkte der Hdhen durch Schm'tt mi't den Sen ten, f_a ist Fu^punkt von h_a, usw: 
#116: pedal points of the altitudes, f_a is pedal point of h_a, . . . : 

#117: solving h_a = a_ for the parameters 
#11S: SOLVE (h_a = a_, [t, tl]) 


#119 : 


# 120 : 


xl- (x2 - x3) - x2 + x2-x3 + (yl - y2)-(y2 - y3) xl- (y2 - y3) + x2- (; 

t = A "tl — 


2 2 2 
x2 - 2-x2-x3 + x3 + (y2 - y3) 


x2 - 2-x2-x3 


f_a := SUB ST 


a_, t, - 


xl- (x2 - x3) - x2 + x2-x3 + (yl - y2)-(y2 - y3) 

2 2 2 

x2 - 2-x2-x3 + x3 + (y2 - y3) 


#121: f_a := 


xl-(x2 - x3) + (y2 - y3)-(x2-(yl - y3) + x3-(y2 - yl)) 

2 2 2 

x2 - 2-x2-x3 + x3 + (y2 - y3) 


2 2 2 

xl- (x2 - x3)-(y2 - y3) + x2 -y3 - x2-x3-(y2 + y3) + x3 -y2 + yl- (y2 - y3) 


2 2 2 
x2 - 2-x2-x3 + x3 + (y2 - y3) 


#122: auf glei'che Wei'se werden die Fuf^punkte f_b und f_c ermittelt / same procedure to obt 
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#122: auf glei'che Wei'se werden die Fuftpunkte f_b und f_c emn'ttelt / same procedure to obtain f_b and f_c 


#123 : f_b := 


xl -x2 - xl-(2-x2-x3 + (yl - y3)-(y3 - y2)) + x3-(x2-x3 + (yl - y2) ■ (yl - y3)) 


xl - 2-xl-x3 + x3 + (yl - y3) 


xl ■ y3 + xl- (x2-(yl - y3) - x3- (yl + y3)) + x2-x3- (y3 - yl) + x3 -yl + y2-(yl - y3) 

2 2 2 
xl - 2-xl-x3 + x3 + (yl - y3) 

2 2 
xl - xl- (x2 + x3) + x2-x3 + (yl - y2)-(yl - y3) xl - xl- (x2 + x3) + x2-x 

#124: f_c := (x2 - xl) + xl, 


xl - 2-xl-x2 + x2 + (yl - y2) 


xl - 2-xl-x2 


Tabelle fur die Fusspunkte der 3 Hbhen / table for the pedal points of the altitudes: 

Fusspunkt der Hbhe h_a: f_a 

fusspunkte := Fusspunkt der Hbhe h_b: f_b 

. Fusspunkt der Hbhe h_c : f_c . 

ped point of altitude h_a: f_a 

pedpts := ped point of altitude h_b: f_b 
. ped point of altitude h_c : f_c . 

Die Schwerli'm'en/Mi ttelli'm'en verbi'nden die Ecken mi't den gegenuberli'egenden Sei tenmi tten , 

ihr Schm'ttpunkt ist der Schwerpunkt/Mi ttelpunkt des Drei'ecks. 

the Medians connect the vertices with the opposite Midpoints of the sides, 

their coMMon intersection point is the centroid of the triangle. 

[s_a := pi + q- (M_a - pi) , s_b := p2 + q- (M_b - p2) , s_c := p3 + q- (m_c - p3)] 

[schwerli'm'en := [s_a, s_b, s_c], Medians := [s_a, s_b, s_c]] 

r 2 2 1 

SOLVE (pi + ql- (M_a - pi) = p2 + q2- (M_b - p2) , [ql, q2]) = I ql = a q2 = 


#136: SUB ST pi + ql- (M_a - pi) , ql, = 


xl + x2 + x3 yl + y2 + 


I" xl + x2 + x3 yl + y2 + y3 

#137: schwerpunkt := , , centroid := schwerpunkt 
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■LOADCD : \DfD\DNL\DNL96\MTH22\TRIANGLE . rath) 


# 2 : 

#3: 

#4: 

# 5 : 

# 6 : 


#7: 

# 8 : 

#9: 

# 10 : 


# 11 : 

# 12 : 

#13: 


#14: 

#15: 


[xl := -6, yl := -3, x2 := 5, y2 := -1, x3 := 1, y3 := 4] 


sen ten 
umkrei'smpkt 
umkrei's 
sei tsymm 



angbi'ss 
incenter 
inci rcle 
[ta, tb, tc] 


hoehen 
orthocenter 
[f_a, f_b, f_c] 


medians 

centroid 


#16: op = 



osculation pts of i nci rcle [x , y] 

on side a a: [2 . 599931251, 2.000085935] 

on side b b: [-0.8109992093, 2.18900079] 

on side c: [1.219986962, -1.687275097] . 

ped point of altitude h_a: [1.682926829, 3.146341463] 

[0.5, 3.5] 

[2.008, -1.844] 


#17: pedpts = ped point of altitude h b : 

. ped poi nt of alti tude h_c : 
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Two Constructions of the 17-Edge 

Josef 

In DNL#20 Johann Wiesenbauer's Titbits dealt with GauB’ proof that a 17-edge can be constructed 
using straightedge and compass only. I’ll show two instructions how to construct the 17-edge. I found 
the first one in Richard Freytag's article 44 Das regulare Siebzehneck ” , DdM 3, 1992. The second one 
goes back to Richmond and was forwarded by Johann Wiesenbauer. The figures were produced with 
WinKon [1], a useful tool for geometric constructions. You can follow the description of the construc- 
tion. 

(krs = circle, ger = line, nor = perpendicular line, pkt = point, mpt = midpoint, str = segment) 

[1] WinKon, Robert P. Michelic, Pillweinstr. 8, A-4020 Linz, Austria 


r=4 

k=krs((0|0),r) 

H(-r/4|0) 

C(0|r) 

c1=str(C,H) 

k2=krs(H,C) 

x=ger((0|0),(5|0)) 

S=pkt(x,k2) 

k3=krs(S,C) 

k4=krs(S’,C) 

T=pkt(x,k3) 

U=pkt(x,k4) 

B(-r|0) 

M=mpt(B,T’) 

k5=krs(M,B) 

F=pkt(k5,ger((0|0),C)) 

r1=|0,F| 

k6=krs(F,r1) 

G=pkt(k6,y) 

k7=krs(G,|0,U , |) 

Q=pkt(k7,x) 

k8=krs(0,|Q,U’|/4) 

A=pkt(k8,x) 

PI 7(r|0) 
h=nor(A’,x) 

P=pkt(k,h) 

s=str(P17,P) 

s2=str(P17,P’) 

l=|P,P17| 

P2=pkt(k,krs(P,l)) 

ns=str(P,P2) 

P3=pkt(k,krs(P2,l)) 

P4=pkt(k,krs(P3,l)) 

P5=pkt(k,krs(P4,l)) 

P6=pkt(k,krs(P5,l)) 

P7=pkt(k,krs(P6,l)) 



P8=pkt(k,krs(P7,l)) 

P9=pkt(k,krs(P8,l)) 

P10=pkt(k,krs(P9,l)) 

P1 1 =pkt(k,krs(P1 0,1)) 

P12=pkt(k,krs(P11,l)) 

P13=pkt(k,krs(P12,l)) 

P14=pkt(k,krs(P13,l)) 

P15=pkt(k,krs(P14,l)) 

P16=pkt(k,krs(P15,l)) 


I'd like to invite you reproducing the constructions using 
any other dynamic geometry program, Josef 
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u_c=krs((0,0),10) 

A(10,0) 

B(0,10) 

D(0,-2.5) 

d_a=str(A,D) 

,<BDA divided by 4 giving (p 
sc=krs(D,|D,A|) 

P1=pkt(sc,y) 

hc1=krs(P1,8) 

hc2=krs(A,8) 

si =ger(D,pkt(hc1 ,hc2)) 

P2=pkt(s1 ,sc) 

hc3=krs(A,5) 

hc4=krs(P2,5) 

hc5=krs(P1,5) 

s2=ger(D,pkt(hc3,hc4)) 

s3=ger(D,pkt(hc5,hc4)) 

a1=ger(D,P3,45°) 

H=pkt(a1,x) 

hc6=krs(mpt(H,A),|A,H|/2) 

K=pkt(y,hc6) 

E=pkt(x,ger(D,P3)) 

hc7=krs(E,|E,K|) 

A1=pkt(hc7,x) 

M=nor(A1’,x) 
l2=nor(A1 ,x) 

.points S3 and S14 
Q=pkt(u_c,l2) 

.points S5 and SI 2 

R=pkt(u_c,l2) 

hc8=krs(Q’,4) 

hc9=krs(R’,4) 

s4=ger(0,pkt(hc8,hc9)) 

S4=pkt(s4,u_c) 

,S4=4 

w1=str(R’,0) 

w2=str(Q’,0) 

d=|R’,S4| 

k1=krs(A,d) 

Z=pkt(k1,u_c) 

Z 1 = p kt( u_c, krs(Z’ , d )) 

Z2=pkt(u_c,krs(Q,d)) 

Z3=pkt(u_c,krs(R,d)) 

Z4=pkt(u_c,krs(Z3,d) 

Z5=pkt(u_c,krs(Z4,d)) 

Z6=pkt(u_c,krs(Z5,d)) 

Z7=pkt(u_c,krs(Z6,d)) 

Z8=pkt(u_c,krs(Z7,d)) 



#1 :A3 

#2:A5 

#3: cp 

#4: cp 

#5: cp 

#6: cp 

#7: cp 

#8: cp 

#9:1 

# 10:2 

#11:3 

#12:4 

#13:5 

#14:6 

#15:7 

#16:8 

#17:9 

#18:10 

#19:11 

#20:12 

#21:13 

#22:14 

#23:15 

#24:16 

#25:17 
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Hello Partitioners 

The discussion about the number of partitions raised by Albert Rich in DNL#20 has provoked 
a fruitful discussion. You can benefit with an improvement of NUMBER. MTH, Josef. 

From: “Soft Warehouse (Albert Rich)” <s wh@aloha.com 20 March 1996 

Hi Josef, 

Thought you might be interested in the following: 

To: Johann Wiesenbauer 
From: Soft Warehouse 

Dear Johann, 
you wrote: 

Attached you will find the file parts, mth containing some implementations of the partition func- 
tion parts(n). I do hope that in particular the last one lives up to your expectations. (If it turns out to 
be inaccurate, adjust all accuracy by increasing 11 in log(n,ll). Iam not sure whether 11 suffices in 
all cases which is certainly a weak point!) Looking forward to your answer. 

Wow! ! ! PARTS( 10099) in 45 seconds is impressive. My hats off to Ramanujan and Rademacher for 
coming up with the formula and to you for implementing it in DERIVE and making it available for 
others to use. I will add it immediately to NUMBER.MTH and give credits to you and R & R. 

Concerning the accuracy required, I will ask Jim FitzSimons, who co-authored the original PARTS 
function with me to compare your PARTS( 10099) with the corresponding result given by Maple V, 
and the time reqired. 

I hesitate to even ask this, but is there a related formula for computing the DISTINCT partitions of a 
number? Perhaps the number of distinct partitions can easily be computed knowing the number of par- 
titions. 

Thanks and Aloha, Albert Rich 

From: “Soft Warehouse (Albert Rich)” <s wh@aloha.com 21 March 1996 

To: Johann Wiesenbauer 

Hello partitioners, 

The following is a printout that I used to check your partitions function: 

PARTS (n) := 

If n < 2 
# 1 : 1 

FLOOR (APPROX (1(1/77 • j (k_/2) - X(I F (GCD(h_ , |<_) = 1, C0S(tt- (X(( i_/k_ - 1/2) - (M0D(i_-h_/IO - 1/2), 


i_, 1, k_ - 1) - 2-n*h_/k_)). 0), h_, 1, k_)-(- 2->/6-e A (- 7i-,/C24-n - l)/(6-k_)) • (e A (7T-,/(24-n - 


l)/(3-k_))-(6-k_ _ 7rV(24.n - l)) - 6-k_ - 7T.,/(24-n - l))/(k_- (24-n - l) A (3/2))), k_, 1, ,/n/L0G(n, 


11)), LOG(l/(4.nV3)-EXP(77-y(2-n/3)), 1Q) + ^ + 0 .5) 
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#2: The following three properties are due to Ramanujan: 

#3: (Hardy and Wright, An Introduction to the Theory of Numbers, page 287): 

#4: M0D(PARTS(5-m +4), S) = 0 

#5: M0D(PARTS(7-m +5), 7) = 0 

#6: MO D( PARTS (11 -m +6), 11) = 0 

#7: We can use them to verify PARTS(n) as follows: 

#8: VECTOR (MOD(PARTS(S-m + 4) , 5), m, 20) 

#9: [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 

#10: VECTOR (M0D(PARTS(7-m + S) , 7), m, 20) 

# 11 : [ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 

#12: VECTOR (M0D(PARTS(11- m + 6) , 11), m, 20) 

#13: [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 

#14: M0 D ( PARTS ( 5 ■ 1000 + 4), S) = 0 


, 0 , 0 , 0 , 0 , 0 ] 

, o, 0, 0, 0, 0] 

, o, o, o, o, 0 ] 


Johann Wiesenbauer told me that MAPLE and MATHEMATICA as well hang up calculating 
parts (10099). One obtains the correct results only after using some calculating time consuming de- 
tours. Josef 

PARTS(n) is now included in the utility file Combi natori al Functi ons . mth of DERIVE 6. 

From: “Soft Warehouse (Albert Rich)” <swh@aloha.com 26 March 1996 

To: Johann Wiesenbauer 

Dear Johann, 

at 01:40 PM 3/25/96 you wrote: 

Concerning the problem of writing an efficient DERIVE-routine for the number of decomposi- 
tions of a natural number into distinct summands I have come across a formula (namely Theorem 354 
in Hardy & Wright's book on number theory) which seems quite suitable to this end. You will find the 
definition together with some examples in the file DPARTS.MTH attached to this e-mail. To make it 
self-contained you should replace p(n) by its definition according to R&R. As you might guess I am 
very eager to learn what you think of it. 

The behaviour of DERIVE 's summation routine changes when the difference between the upper and 
lower limits is greater than 25 because above 25 DERIVE attempts to compute an antidifference for 
the sum. If that succeeds, the antidifference is used to compute the sum directly instead of actually 
summing up the terms. This can be a big win if a large number of terms are to be computed and/or it is 
difficult to compute these terms. Note that as of DERIVE and DERIVE XM version 3.1 1, the differ- 
ence in the limits before the antidifference algorithm was used raised from 25 to 100. 

Unfortunately, as you discovered, computing the antidifference an lead to problems. The summand 
must be simplified before the antidifference can be computed. In your definition for PD0, you use the 
predicate 


NUMBER (SQRT (8*n-16*i +1)), 
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in the IF construct of the summand. This predicate simplifies to “false” since i_ is not known at the 
time the summand is simplified and the resulting expression is not a number. Therefore, the IF con- 
struct returns the else clause argument, which in this case is 0. 

You avoid this problem by making a vector and the summing its elements. To avoid making and hav- 
ing to store the vector, the above predicate can be replaced by 

8*n-l 6*i_+l = FLOOR(SQRT ( 8*n-l 6*i_+l ) ) A 2 

which does not simplify to “true” or “false”. Therefore, unless you see any problems, I propose the 
following definition for NUMBER . MTH. 

DISTINCT_PARTS (n) := 

SUM (IF ( 8*n-l 6*i_+l=FLOOR (SQRT ( 8*n-l 6*i_+l ) ) A 2, PARTS ( i_) , 0 ) , i_, 0 , n/2 ) 

The bug in my original version of DISTINCT_PARTS (and PARTS for that matter) is more subtle and 
difficult to revolve. It arises because a variable name conflict occurs when the summand in a user- 
defined function recursively calls a function. I know of no easy way to avoid this problem. 

Thus I greatly appreciate your partition functions not only because they are fast but because they avoid 
this variable name conflict problem. 

Aloha, Albert Rich, Applied Logician 


From: “Soft Warehouse (Albert Rich)” <s wh@aloha.com 1 May 1996 

To: Johann Wiesenbauer 

Dear Johann, 

As always, I was amazed by your DERIVE programming prowess shown in your Titbits column in the 
DNL#21. 

I made some changes to your iterative definition of LUCAS that takes advantage of some features of 
DERIVE that simplify and clarify the definition somewhat. The following is my definition: 

LUCAS (n):= ( ITERATE ( IF ( (n AND d_) =0 

[a_ A 2-2* (-1) A c_, a_*b_- (-1) A c_, 2*c_, d_/2 ] , 

[a_*b_- (-1) A c_,b_ A 2+2* (-1) A c_, 2*c_+l , d_/2 ] ) , [a_,b_,c_,d_] , 
[2,1,0, 2 A FLOOR (LOG (n, 2 ) ) ] , FLOOR (LOG (n, 2 ) ) + 1)) SUB 1 

Note the use of the AND operator to perform the bit-wise logical AND operation on integers (see Sec- 
tion 4.16). This makes it possible to use only a four-element vector for the iteration. 

Also 2 A FL00R(L0G(n , 2) ) is a more direct way to compute the largest power of 2 less or equal to n. 
Aloha, Albert Rich 


From: “Soft Warehouse (Albert Rich)” <s wh@aloha.com 4 May 1996 

To: Johann Wiesenbauer 

Dear Johann, 

at 08:16 PM 5/4/96 you wrote: 

Thank you very much for your words of praise as regards my programming efforts (it is a good 
feeling to know that there are people out there like you who care for them) as well as your valuable 
improvement of the function LUCAS(n), which certainly looks much tidier now. If you don 7 object to 
it I will publish it in my next ‘ Titbits 
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Please do! It is a good way to highlight DERIVE Version 3’s ability to do bitwise logical operations. 

Thank you also for your e-mail next to the last, although it was a bit shocking for people like me who 
used to think that DERIVE was relatively bugfree. 

Dave and I have been “teaching” mathematics to muMATH and DERIVE for the last 18 years. In that 
time we have built up a test suite of over 9000 problems that DERIVE must pass before a new version 
is released. However, bugs will always exist in large, complex systems like DERIVE that attempt to 
automatically simplify any mathematical expression thrown at it. 

Therefore, one should never trust the results of any computer algebra system, including DERIVE. We 
highly recommend verifying results produced by DERIVE (e.g. substitute solutions of equations back 
into the original equations; differentiate antiderivatives, etc.). 

Well, of course I have been to used encountering bugs every now and then like e.g. 

LIM(x • EXP(- x a 2) • INT(EXP(t A 2) , t, 0, x) , x,inf) = 0 

(cp. LIM(DIF(x • INT(EXP(t A 2) , t, 0, x) , x)/DIF(EXP(x A 2) , x) , x, oo) = 1/2) 

Thanks for pointing out the above bug. I have fixed the problem and will send an update to Soft 
Warehouse Europe. Please report any other bugs in DERIVE that you encounter so they can be fixed. 

...but it is much harder to put with the possibility that such common functions like SUM (and supp- 
sedly also PRODUCT) might go haywire out of the blue after lulling the user inio a false sense of se- 
curity by yielding a certain amount of correct values! Maybe there is some remedy for this serious 
problem. 

Using antidifferences (like antiderivatives) is always dangerous, but I do not know of any alternative. 
Perhaps we could have an option to force SUM to not use the antidifference if the lower and upper 
limits are known. 


The Lorenz attractor produced with DERIVE, converted into an Acrospin file and then ani- 
mated. (I will show how to do this in Bonn, Josef) 
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This is a second letter from Alfonso J. Poblacion Saez. He wrote: 


Dear Josef, 

I send you an article which I had thought and written almost a year ago. I did not send it to you be- 
fore, because in fact it doesn't deal with anything especially relevant from the mathematical point of 
view; it is only a curiosity. Finally, I dare to send it to you. 

Another question related: I was thinking (a very dangerous thing, you will see) why do not have a 
brief section in the newsletter in which DUG members can propose puzzles, problems (even unsolved) 
or mathematical recreations for which DERIVE can be applied to solve them. Some contributions pub- 
lished fit perfectly this idea like Conway's Game of Life, some beautiful drawings, the tennis net 
analysis, etc. It could be named, for example as AC DC (you know, the Amazing (or Amusing) Comer 
of the DERIVERs' Curiosity), or any similar thing. 

Forgive me if I am completely out of order, it is only a suggestion. 

With the four pages of the squaring circle article I enclose some DERIVE sheets. They are not to be 
included. It is only to show you how I did. Perhaps they will have better room in the diskette of the 
year. 

Yours faithfully 
Alfonso J.P.S. 

This is a great idea. As you can see, we have the AC DC section. We will start with your ap- 
proximations of squaring the circle. I could imagine that students will like your paper and will 
be eager to find the solutions. Josef 

Some Approximazions to Square the Circle 

Alfonso J. Poblacion Saez, Valladolid, Spain 


INTRODUCTION 

Recreational Mathematics are generally considered only as an entertainment to kill free time. 
However, sometimes a lot of puzzles and challenging problems become interesting and it is necessary 
to apply relevant results to solve them. 

One of the most popular problems since the Greeks has been the squaring of a given circle. In 
1882 Ferdinand Lindemann proved by algebraical methods that this question has no solution because 
n is transcendent and hence is not constructible. But a great deal of people lent us some beautiful and 
ingenious approximations. In this brief article I have gathered four of them (from [2], [4] and [6]) in 
order to find ist accuracy using DERIVE. Although it is only a recreation, perhaps someone could use 
them as easy exercises in analytic geometry or to introduce some aspects of a part of the mathematics 
usually forgotten: the History of Mathematics and the evolution of the problems (see [1] and [3] for 
further details in this sense). So, if you have some spare time, try to answer the questions proposed! 

THE PROBLEM 

To square a circle consists in constructing a square having the same area as a given arbitrary cir- 
cle, using only straightedge and compass. This last condition means that we can only do these things: 
to choose points in the plane, to draw a straight line passing through two points, to draw circles of a 
given center and radius and to find intersection points of lines, circumferences and lines and circum- 
ferences. 
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HOW TO USE DERIVE 


Use it as you like. I have defined some functions to automate the drawings and to find the alge- 
braic resulats, for example: 


SEGM(x,y) := 

MID(x , y) : = 
EC_LINE(x_,y_) := 
DIST(x,y) := 


to plot segments [x, y\ in connected mode, 
to give the midpoint of a segment [x,y], 
to find the equation of a line passing points x_ and y_, 
to calculate the distance between two points x and y. 


and so on. For simplifying the constructions you can begin with a square and a circle with sides of 
length one unit, respectively. 


1. One of the oldest approximations 

The following procedure gives one of the best rational approximations to n and also one of the 
oldest. It was first used by a Chinese astronomer Tsu Ch'ung-Chih who lived in the 5 th century (pic- 
ture on page 50). 


Let ABCD be a given square (side length 1). Take E 


such that AE = — AD and then F on BE verifying 


BF = 


AB 


Consider the perpendicular to AB through F, giving 
point G. Finaly find H such that the segment FH is 
parallel to EG. 


7i can be approximated by 3 + HB. 

What is the very known rational value and how many 
digits are equal to 7 i? 



2. How wrong was Hobbes? 

Thomas Hobbes (1588-1679) (picture on page 50) was mainly a philosopher (remember his 
4 homo homini lupus') but when he was 67 he took an interest about mathematics. He criticized those 
who applied algebraic methods to geometry, like John Wallis. Both of them were involved in a hard 
and long dispute. Hobbes found a dozen methods to square the circle. Here is one of the best in which 
you can use analytic geometry, in spite of Hobbes. 

Starting again with a given square, draw circular 
arcs AC (center D ) and BD (center A). F is such 
that arc CF = 0.5 CE. Then take the distance 
from F to segment CD and obtain G so that FG is 
equal to that distance and parallel to side BC. 

Draw line EG and find H = EG n BC. 

Hobbes assured that the circumference contains 
the arc CE twelve times and that CE = CH , and 
so 7i is approximately 6-CH. 



How wrong was Hobbes? 
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Clear Function Parameters Representation 

Sergey V. Biryukov, Moscow, Russia 


This report is about clear function parameters representation & ode_appr.mth extension. 

EULERO function which has only 1 st order accuracy was included in ode_appr .mth for educational 
purposes. I think that some purposes require this function extension to the system of ODEs and 
modification to the 2 nd order accuracy. Appropriate functions are defined in lines #4 and #5. 

Functions with a numerous and long vector of parameters are difficult to percept. RK( r , v , vO , h , n) , 
for instance which solves the system of ordinary differential equations has three first vector argu- 
ments and two last scalar ones. It would be better to arrange all data in one or several matrices and 
include labels for rows and/or columns. It can be done by Interface Functions which convert clear 
matrix-form arguments to the list of arguments, call the active function and return the result with 
appropriate labels. Such functions can greatly improve clearness of DERIVEed expressions and thus 
will be helpful for teaching. So, Interface Functions Style must be established. This report is the first 
attempt. Expressions #17 and #18 seem to be more alike mathematical notation than #16. Further- 
more they are commented and thus more suitable for teaching and using by novices in DERIVE. On 
the other hand #16 is more compact and clear for design and analysis by professionals. 

Methods which give different orders of accuracy are compared by the example of one dimensional 
motion (#12 - #15, #20 - #25). 


#1: "=== ODE_APPR . MTH Extension for ODEs Systems === ODEAPPR1 . MTH === 7 Aug 94 ===" 

#2: "Uses ODE_APPR . MTH & VECTOR. MTH Copyright (c) 1994 by Sergey V. Biryukov" 

#3: [E (v, i) :=ELEMENT (v, i) , EE (v, i, j ) :=ELEMENT (v, i, j ) , D (v) : = DIMENSION (v) , 

D_E (v, i) : = DELETE_ELEMENT (v, i) ] 

#4: EULER_SYS (r,v,v0,h,n) : = ITERATES (v_+h • APPEND ( [ 1 ] , LIM (r , v, vj ) , v_, vO, n) 

#5: EULER_2_SYS (r,v,v0,h,n) : = ITERATES (v_+h • APPEND ( [1] , LIM (r, v, v_+h/2 • APPEND ( [1] , 

LIM (r , v, v_) ) ) ) , v_, vO , n) 

#6: " ===== Clear System of ODEs Representation (Notation) ===== " 

#7 : EULER_SYS_(m,p) :=APPEND ( [VECTOR(LHS (EE (m, i , 2 ) ) , i , 2 , D (m) ) ] , 

EULER_SYS (VECTOR (RHS (EE (m, i , 1 ) ) , i , 3 , D (m) ) , VECTOR (LHS (EE (m, i , 2 ) ) , i , 2 , D (m) ) , 
VECTOR (RHS (EE (m, i , 2 ) ) ,i,2, D (m) ) , EE (p, 2 , 1 ) , EE (p, 2 , 2 ) ) ) 

#8 : EULER_2_SYS_ (m, p) :=APPEND ( [VECTOR(LHS (EE (m, i , 2 ) ) , i , 2 , D (m) ) ] , 

EULER_2_SYS (VECTOR (RHS (EE (m, i , 1 ) ) , i , 3 , D (m) ) , VECTOR (LHS (EE (m, i , 2 ) ) , i , 2 , D (m) ) , 
VECTOR (RHS (EE (m, i , 2 ) ) , i , 2 , D (m) ) , EE (p, 2 , 1 ) , EE (p, 2 , 2 ) ) ) 

#9 : RK_ (m, p) :=APPEND ( [VECTOR (LHS (EE (m, i , 2 ) ) , i , 2 , D (m) ) ] , 

RK (VECTOR (RHS (EE (m, i , 1 ) ) , i , 3 , D (m) ) , VECTOR (LHS (EE (m, i , 2 ) ) , i , 2 , D (m) ) , 

VECTOR (RHS (EE (m, i, 2 ) ) , i, 2 , D (m) ) , EE (p, 2 , 1) , EE (p, 2 , 2 ) ) ) 

#10 : [E2C_(m, i, j ) : = D_E (EXTRACT_2_COLUMNS (m, i, j ) , 1) , 

"Delete 1st Row & Extract 2 Columns"] 


# 11 : 
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The CHAOS GAME on the TI-92 

(by Larry Gilligan a.o.) 

In the last DNL I put L.Gilligan's, J. Rose's and N. Rich's book "Mastering the TI-92" on our 
book shelf. In the meanwhile I had the occasion to work with this book and I can highly rec- 
ommend it to you all. You will find a lot of tricks and examples how to use the TI-92. Do you 
know how to create you own tool bars with customized pull down menus? I know because 
I’ve learned it from "Mastering". I want to show one example how to program with the Tl. I 
had not time to ask Larrie & Co for permission but I hope they will understand my idea to do 
some promotion for their product. 

I will demnstrate how to program the well known CHAOS GAME using the TI-92. I changed 
the program from Mastering by adding a short dialogue for the input of the parameters which 
are constant in the original version (the coordinates of the starting point). If you don’t know 
the CHAOS GAME then be patient please, in one of the next DNL's you will find a description 
in Maria Koth's contribution "Computer Graphics with DERIVE". 

Many thanks to Larry & Co in advance. 

The program: 

s 1erp( ) 

Prgm 

Local n.fx.fy, 1 ,x,y , v 
ClrlO 

Input "Initial x-value -6<x<6",fx 
Input "Initial y-value 0<y<4”,fy 
Input "Number of Iterations” ,n 
Cl rGraph 
Cl rDraw 
FnOff 

enext 2 lines set graph window 
-6.5-»xm1n: 6.5+xmax: l-»xscl 
-0.5-*ym1n: 4.5->ymax: l-»yscl : 2->xres 
setGrapht "Axes" , "Off" ) 
f x-»x :© x-value of Initial point 
fy-»y:e y-value of Initial point 
For 1,1, n 
rand(3)->v 
If v=l Then 
( x+0) /2-»x 

(y+4) /2-»y 
Elself v=2 Then 
( x+6) /2-»x 
(y+0) /2-»y 
Elself v=3 Then 
( x-6) /2-»x 
(y+0)/2->y 
Endlf 
PtOn x,y 
EndFor 
EndPrgm 




MAIN 


RAD EXACT 


FUNC 








Finding Integral Curves with the TI-92 

(by Josef Bohm) 

In the last weeks I had to teach Calculus in one of my classes. It is unevitable to solve the 
problem of finding an integral curve passing a given point. This was a great occasion to 
demonstrate the Tl's abilities. 

2x 

Example"!: Let's find the integral curve to /' = e 5 x + 1 passing P(2|-3). 

You can follow the screen shots. There is a little trick which is not visible. If 
you edit /(f,x,c) then you will obtain the antiderivative of fwith respect to x and 
an integration constant c. 


hi'-Twdf t r-t T > ■! \ 

II - f~- h : ^Lsr.L-lz. f * t Prsn 1 0 CTI ■■ 




» &clue 


■■/■Tdx 



“2“ 

4 + X + C 


■ WUU 

:■ -3 = /f<!x 1 

x s 2 


•3 = c--S-L +2 


5-x-e 

9 


lc<f,x,2, ~3> 


m:m .... nm emao func ssiq 




The more advanced user will try to calculate a function i c which does his work. We open the 
Program Editor and edit a new function called i c: 


fFiT^Y F2» ' 

'•'f— (Control 

FE Y Fh’T Y ”1, 

I/O Uar Find. . . Mode 

sic(u,u,x0,yi 
: Func 

:TCu, w)+y0-(, 

: EndFunc 

35 ’ ” 

fCu, V) 1 O=X0> 


And the function does the expected work. 
But there is one more way to solve the 
problem: use Solve from the F2-Menu and 
build the full construction in order to solve 
the equation for the unknown parameter c. 

Example 2: y' = x 2 ; P(3|4) 

Example 3: ^ = t-e~ at ; P(4|4) 


j i yk br * \ c ' |Frgn rolC l % *:■ 


■ ic(x 2 , x, 3 , 4 ) 

■ = J(x 2 )dx | x = 3 arid y = 4 3 c) 

c = "5 

j* solve(y = a ^t)d J L|t = 4 and y = 4 > c) 
_ e ^ a + 4 

ic<e A < ~a*t >*t .t.4,4) 

MAIN ftflD *UTP FUHC Fr'IB 


- £ K1 j ! C |gj -L ^|P * g n 1 0 P =: ■■ 

c - “5 

* = f (sf a ^ ] i = 4 and y=4*c] 


TTwjy rjz v r : T T ** T 1;i 

T t ; ■ Pr'-^rilO 


-JO ^ 't)dt 



lic<re .4,4) 


IHMH ftliTII FUNC &2U 


r pii: w ni H n u tv? func i th 
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Classic DERIVE compared with TI-92 - DERIVE 

(by Josef Bohm) 


The Financial Mathematics Superformula 


Long, long ago, in 1991, I wrote an article on financial mathematics (DNL#1,#2, 1991). I dis- 
cussed the use of the financial functions offered by DERIVE v.1. In the TI-92 manual I found 
an application (App. 15) which reminded me on these functions and gave the impact to im- 
prove the function presented there to find payments, future & present net values, the number 
of payments and the interest rate for the TI-92 and to use also DERIVE v.3's new capabilities 
to build a general financial mathematical function for all purposes. 


See first the classic DERIVE version: 

(You can find the respective text for the examples on the next page.) 

#1: [Precision Approximate, Notation := Decimal , PrecisionDigits 15] 

#2 : InputMode : = Word 

#3: H (i, ip, rp) := (1 + i/ (ip • 100) ) A (ip/rp) 

#4: FMF (cap, f l_p0 , paym, n, el_b0 , i , ip, rp) : = FLOOR (100 -RHS (NSOLVE 

(cap=paym- (H ( i , ip, rp) A n-1 ) / (H ( i , ip, rp) -1 ) • (H (i, ip, rp) A 
( l-el_b0 ) /H (i, ip, rp) A (n- (l-fl_p0) )),x,0,~))+0.5)/100 

"The examples:" 

FMF (15000, 0,x, 48,1,9, 12 , 12 ) =373 . 28 
FMF (x, 1,100, 20, 0,5. 5, 2, 4) =2 313. 19 
FMF (100000, 1, 1500, x, 1, 6, 2, 2) =37 .17 
FMF ( 114720, 0, 5439, 24, l,x, 1,12) =13. 4 9 
FMF (x, 0,2500, 15, 1,3. 75, 1,1) =28288. 24 
FMF (x, 1,1500, 37, 1,6, 2, 2) =992 61. 33 

And now let us try with the TI-92!: 

We open the Program Editor and then we edit the following program: 

: fmf ( tl ,sl,t2,t3,s2,t4,t5,t6) 

:Func 

: Local tl ,s1 , tempfunc, temps trl ,auxl 
: (l+t4/(100*t5) ) A (t5/t6)-»auxl 

: -tl+t2*(auxl A t3-l )*auxl A (l-s2) / ( (auxl-1 )*auxl A (t3-t3*sl) )-»tempf unc 
-.For tl ,1 ,4,1 

: H t"&exact(str1ng(t1) )-»tempstrl 
: I f when(#tempstrl=0 , false, false, true) Then 
: If t 1=4 

:Return f 1 oor ( 100*approx( nSol ve( tempf unc=0 ,#tempstrl ) | #tempstrl>0 and 
#tempstrl<100)+.5)/100 

: Return f 1 oor ( 100*approx( nSol ve( tempf unc=0 ,#tempstrl ) ) + . 5) / 1 00 
: Endlf 
: EndFor 

:Return "parameter error" 

: EndFunc 
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Explanation of the variables used: 

cap = tl amout of money (present value or future value) 

f1_p0 = si 1 if cap/tl is a future value or 0 if cap/tl is a present value 

paym = t2 periodical payment(s) 

n = t3 number of payments 

e1_b0 = s2 1 if the payments are due at the end of the periods 
or 0 if they are duw at the begin of the periods 
i = t4 annual interest rate (nominal) 

ip = s3 number of interest periods/year 

rp = s4 number of payments/year 

The Examples: 

1 . Find the monthly payments for a loan of 


3. A deposit of ATS 100000 should be gathered by payments of ATS 1500 twice a year at the end of 
the semesters. The interest rate is 3% per semester. How many full payments are necessary? 

4. A car is priced at ATS 144 720. You have to pay an account of ATS 30 000 and the remaining 
balance in 24 monthly payments of ATS 5439 (at the and of each month). What is the annual in- 
terest rate? 

5. Similar to example 2. 

6. Which amount of money is necessary that you can spend 15 annuities of 2500 each at the end of 
the years if you gather 3.75% per year? 


$15,000 at an annually interest rate of 9% 
with twelve interest periods/year. The loan 
is paid back within 4 fours at the end of the 
months. 

If I would save quarterly at the beginning of 
each period $100 through 5 years at 2.75% 
per semester, how much money would I end 
up with? 


^ HfilqebralCalc iQtfoer |FrgnIQ|Clear a-z... | 

1 f nf C 15000 , D , x ? 48 , 1 ? 9 , 12 , 12) 373. 2§ 

ifmftx, 1 , 100 ? 20, 0, 5.5, 2, 4) 2313.19 

i f IO00OQ , 1 , 1500 , x, 1 , 6, 2, 2) 37. 17 

1 frcf( 144720 - 30000 ? 0 , 5439 i 24,1,k,M> 

13.49 


■ f nfCx , 1 , 1500 , 37 , 1 , 6 , 2 , 2) 

2500 j 15, 1, 3.75, 1, 1) 

fnf<l 

MHIN KM A urn FUHC 9/30 


99261.33 


28288.2 
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Jose Verhoosel 

Using DERIVE 3.111 encountered the next problem: 

Using DERIVE 2.5 1 it is possible to solve a • [[1,2], [3, 4] ] = [ [1,1], [2,3]]. The result is the matrix 
a, which fits the equation. This is not possible any more within DERIVE 3.11. Can you give me the reason or is 
this just another bug? 

DNL: Albert Rich 's answer is: Beginning with DERIVE 3, variables can be declared integer, real, complex or 
nonscalar. By default variables are assumed real. 


Josef Lechner, Viehdorf, Austria 

Has anybody of you ever tried to find the general solution of a quartic using DERIVE ? 
sol ve(ax A 4+bx A 3+cx A 2+dx+e=0 , x) 


Modem didactics and Pedagogy? 


“Teach your scholar to observe ... you will soon raise his curiosity. Put the 
problems before him and let him solve them himself. Let him know nothing 
because you have told him, but because he has learned it for himself. Un- 
doubtely the notions of things thus acquired for oneself are clearer and 
much more convincing than those acqired from the teaching by others ... 


No, this was written by Jean-Jacques Rousseau nearly 260 years ago: Emile (1762) 




Tsu Ch'ung-Chih (429 - 500) 


Thomas Hobbes (1588 - 1679) 


http://hua.umf.maine.edu/China/astronomv/tianpage/0014ZuChongzhi9296bw.html 


http://oregonstate.edu/instruct/phl302/philosophers/hobbes.html 




